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ARSTRACT 


Carmichael fijJ has proposed a direct injection, 
temperature ignition (T1), reciprocating engine design which 
divides the conventional four stroke cycle functions between 
a compression cylinder and an expansion cylinder, 
imterconnected by a static regenerative heat exchanger. The 
prediction of heat transfer loss in the expansian cylinder is 
required to predict the performance of the new cycle. 


Woschni’s correlations [9], developed for direct 
injection, quiescent, four stroke, compression ignition (CTI) 
engines, are used to predict uncorrected expansion cylinder 
neat transfer loss. laws of similarity allow, in principle, 
extrapolation outside the experimental data range for forced 
convection heat transfer only. General applicability of 
Woschni ‘’s correlations is established by e@xamining the 

xpansioan cylinder quiescent condition, using predicted mean 
gas velocity magnitudes and fully developed, turbulent flow, 
forced convection models. Fredicted radiant peak flame 
emissivity and temperature are found and compared with 
nominal peak values observed in conventional CI engine cycles 
to establish a minimum radiation flux reduction ratio. I" 
mean reduction ratio for forced convection flux is alsa 
found. A corrected, conservative prediction of heat transfer 
loss in the expansion cylinder is then established. 
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CHAF TER ONE 
INTRODUCTION 


i SeOrE 

Conventional compression ignition ee reciprocating 
engines have found significant application in commercial, 
industrial, and in land and marine transportation systems. 
Recent emphasis on fuel costs and avVallability, in 
conjunction with improved materials engineering, has prompted 
consideration of design alternatives which offer increased 
thermal efficiency and reduction in operating costs. A New 
Cycle Engine design has been proposed by Carmichael Cid in 
which the functions of a conventional four stroke 
reciprocating cylinder are divided between a compression 
cylinder for charge induction and compression, an expansian 
cylinder for charge combustion-expansion and exhaust, 
interconnected by a static regenerative heat exchanger. 

The attractiveness of such an engine 1s possible 
because improved fabrication techniques presently allow 
construction of strong, thermal resistant, extremely thin- 
walled ceramic matrices possessing high heat transfer 
properties. A cyclic interchange of compressed charge into 
and expired exhaust gas out of the expansion cylinder via the 
regenerator establishes an effective medium for heat transfer 
and allows use of exhaust heat otherwise lost to the system. 
The temperature of charge air heated By adiabatic compression 


is Boosted from the regenerator heat exchange, a direct 





consequence of which 1s significant reduction ian cylinder 


required to achieve optimum 


compression mMatio, es 


efficiency. A computer-generated thermodynamic model of this 
engine in preliminary design indicates that improved thermal 
efficiency may be achieved. An important element in the 
cycle analysis is accurate estimation of heat transferred 
from the gas medium in the expansion cylinder to its 
enclosing surface area, and ultimately to the coolant. The 
magnitude of this heat transfer reduces both heat 


availability to the regenerator and the thermal etficiency. 
The purpose of this thesis 1s as follows: 


We To evaluate the applicability of existing 
heat transfer correlations developed to 
predict total mean instantaneous heat 


transfer coefficient, Hh in four stroke 


ea) 
direct injection, quiescent, CI engines 
for prediction of total heat transfer in 
the New Cycle Engine expansion cylinder. 

Ze To evaluate and compare the magnitude of 
the reduction ratio in both the forced 
convective and the radiant heat transfer 
from the high peak combustion pressure 


conditions inherent in high ic et 


Gon 


engines to low peak combustion pressure 


conditions inherent in the low Fr New 


c 7 


Cycle Engine design. 


ay 





if 


: To use the reduction ratios to correct 
existing heat transfer correlations, at 
applicable, and establish a correction 
allowing better prediction esi heat 
transfer loss in the New Cycle Engine 
expansion cylinder. 
ee NEW CYCLE ENGINE DESIGN BACKGROUND 
A vertical cross-section of the New Cycle Engine design 
is provided in Figure i-l. For simplicity, each cylinder is 
considered of simple geometry, with straight: forward walls 
and flat cylinder head/piston crown configuraticans; bore 
diameter, 8B, equals stroke length, L. Gas velocities are 
boosted by piston speed, imeateGuion SUGELOnN, Ccharge impulse 
ec cme Fegenerator port, .and combustion turbulence only: 
quiescent flow should predominate with no deliberate design 
feature incorporated to accelerate any rotational velocity 
incurred by off-set regenerator port impulse. The cylinder 
crank shafts will be mechanically linked to enable tne 
compression piston to cycle at twice the speed of the 
expansion Histon, alternating between an lmelictivon— 
compression and a free motoring stroke. Proper crank angle 
phase shift will exist to sequence the compression piston 
motion at the start of expansion piston down-stroke to 
optimize charge transfer. 
Correct value timing 1S imperative and is seqguenced 


through one charge cycle as foOilows: 





1 Compression cylinder induction stroke: 
induction inlet Valve 1 Opens; 
regenerator transfer valve <& is 
closed. 

0 Compression cylinder compression 
stroke: Valve 1 closes; Valve <& 
remains closed. 

aie Expansion cylinder combustion- 

xpansion stroke: regenerator exhaust 
Value 2 closes; valve 2 opens inducing 
compressed charge impulse through the 


regenerator as the expansion piston 


e starts the down-stroke; valve 2 then 
closes. 

4. Expansion = vaeedee “Haust stroke: 

valve 2 opens, regenerating heat 

xchange capacity during exhaust gas 

x#piration. Valve 2 remains closed; 


Valve 1 is cycling to support the next 
charge cycle as expansion piston 


reaches top-dead-center (TDC). 


A vertical cross-section schematic of the xpansion 


cylinder and piston linkage showing principle dimensions is 


provided im Figure 1-2, A derivation oa geometric 
relationships required to calculate instantaneous xpansioan 
cylinder heat transfer surface area, Aut volume, VV, and 
piston speed, Sp, is enclosed as Appendix &. Instantaneous 





expansion cylinder pressure, F, mean gas temperature, Tg, and 
volume, V, at S° crank angle increments from the computer- 
generated thermodynamic model is incloséd as Appendix A, 
Table A-i Cid: associated F-V and Ta-V diagrams for the 
entire cycle are illustrated in Figures 1-2 and 1-4 
respectively. Other principle desiaqn and operating 
parameters for the expansion cylinder required to support the 
Neat transfer calculations are listed in Appendix A. 
Appendices A and & provide the data point basis for the heat 
transfer calculations. A consideration of piston ring land 
area, cylinder sleeve clearances, and clearance volume, Vc, 
at TDC is neglected. All the area and volumetric 
calculations ae based on simple cylinder geometry. Within 
the accuracy of reciprocating engine heat transfer 
correlations in general, errors associated with these 


assumptions are considered to be insignificant. 


=O 
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CHAF TER TWO 


FORCED CONVECTION HEAT TRANSFER IN THE 


NEW CYCLE ENGINE EXPANSION CYLINDER 


Ee THEORETICAL BACKGROUND 

When a temperature difference exists petween aie solid 
surface and a moving fluid, thermal energy 1s transferred By 
the joint process of conduction and convection. Shear 
stresses acting on the fluid will slow the fluid immediately 
adjacent to the surface, forming a thin boundary tlayer of 
thickness 5. The heat transferred between the bulk fluid and 
the body is transported by conduction in the boundary layer 
as a function of fluid movement source (natural or forced), 


flow character (laminar or turbulent), physical properties of 


the fluid, and body geometry and roughness. The heat 
transfer rate, ce by conduction in any direction x at any 
instant is determined using Fourier’s Law By a heat 


conduction coefficient, EF, and the temperature gradient: 


Q = -KAIT ee 
ox 

In laminar flow, fluid particles flowing parallel to 
the body surface conduct heat in a normal direction. At a 
certain transition velocity, turbulence is created where 
particles Rave a velocity component normal to the Bulk stream 
direction; the intense mixing of individual layers promotes 
rapid transportation of heat to/from the surface. Even ain 


Curbulent flow, however, a thin boundary layer exists due to 








fluid viscosity across which laminar flow and heat transter 
Characteristics are maintained. In turbulent systems, the 
heat transferred between a body surface and moving fluid 
remains limited by Boundary layer thickness, §6, and fluid 
thermal conductivity, FE. 

An ideal linear velocity and temperature gradient 
across the boundary layer between a gas medium and a boady 
surface allows approximation of the temperature gradient as 
CTG-TW) 7/88 substitution into equation (1) yields the 
Newtonian expression for steady flow conditions: 


QO = -KAvT 
Ox 


-KACTQ—-Twl 
¢ 
d 


—-KCTQG—-Tw) 
¢. 


= QO 


where letting h = -K 
¢. 
oO 


Q = hA(Tg-Tw) (2B) 


For real boundary layers of varying velocity and temperature 
Gradients, hh may be determined in some cases by integration 
of the differential equations describing the heat transfer 
process. For more general and practical applications, 
however, laws of Similarity may be used to predict the value 
of fh. According to the laws of similarity, it is sufficient 
to determine the characteristic relations between independent 
dimensionless groups obtained, recognizing that any such 


relations will apply to all similar heat transfer processes. 


aaeat 
a 
ae 





From the Buckingham pi theorem which allows reduction of 
initial variables to minimum dimensionless groups, the 
following dimensionless numbers are used to describe forced 


convection heat transfer ina fluid medium to a solid body: 


ie Nusselt #, Nu = Hl 7k (3) 
ae Reynolds #, Re = Ul _/y, (4) 
ate Frandtl #, Fr = Vi Se = po 7K. (>) 
wheres Vi oe W/p = MOMentum = diffusivity 
(kinematic viscosity) (6) 
« = B/C Pp = thermal diffusivity (7 
M = dynamic viscosity 


p = density 


0} 
I 


specific heat at constant pressure 
l_ = Re characteristic length 
and in general form: 


Ru = function of Re, Fr = ff (Re,Fr) (3) 





In fully developed turbulent flow systems, analysis of heat 
transfer processes 15 enhanced by a simple analogy between heat 
jand momentum transfer which allows prediction of Hh from knowledge 
lof shear stress (friction) acting on the fluid adjacent to the 
wall. Enown as the Reynolds-Colburn Analogy, it is described by 
Rohsenow in the following manner (€2,p.1851: 

"the heat transferred to a surface divided by the 
maximum amount which could be transferred in 
Bringing the fluid to the temperature of the 
surface equals the momentum transferred in 
passing over a Surface divided by the maximum 
amount of momentum which can be transferred in 
stopping the fluid relative to the surface." 


This analogy, exactly true in parallel flow when Fr = 1) and 


reasonably accurate at Fr values not very high or very low, 


7 





may be used in either constant Heat flux or wall temperature 


conditions and may be expressed mathematically as: 


LS 


Mee oletenet acter = Str — = 72 (F) 
where: Stanton #, St = H/p Cou CG) 
oes t/(pu-72) = friction coefficient Cys) 


t apparent shear stress acting on the fluid 


at the wall 


Significant flow experimentation has determined that for low 
turbulent flow conditions over a smooth plate at a distance x 


from the leading edge: 


E., .0592 Re "= (12) 
flat plate “A 


Similarly, for turbulent flow conditions in a cylindrical 


pipe of diameter Ds: 


oa naa Rep °* (15) 
pipe flow 


Substitution of these relations into equation (9) allows for 


the determination of Rh: 


1? flat plate of length 1: 


fe 
fee = = ,O5392 (Re 7*=“)/2 
ee 
PP 
eaten ce U(Re,)~** Pr 2/5 by 
2 - - 
ON 





substituting for density, the equation (5) identity, p = W/V S 


: , ete co = fi 
h, = .0296 K Ux wC, (Re? Fp 
: ais 
recognizing equations (4) and (3) identities Re, = Ux and 
Vi. 
Pr = wC /K s 
rr AC, 
h, = .0296 K Re *? pri’ 
integrating over length 1 to determine mean h, : 
hy =i (i 292976 k Cees pr i’s di 
1 de x 
0296 K (UZy, 1°" Pro? |) xe * ax 
t ie 
eee kk (Use Prd iron 
os 1 ) : 
Be = .OS70K(Re,)°° Pri (14) 
1 


A similar substitution of equation (13) for turbulent flow in 


a cylindrical pipe yields: 


Np» pipe flow of diameter D: 


a, = -O23 Ethap)*® as (15) 


These particular relations form the basis Gf many heat 
transfer correlations commonly used, rigorously applicable in 
fully developed, steady turbulent flow within estadlished Re 
and Fr limits where temperature difference, Al, across the 


boundary layer 15 not extreme. 





Peete LING CONVECTIVE HEAT TRANSFER IN RECIFROCATING ENGINES 


In reciprocating engines, the significant process by 
which termal energy is lost is By forced convective heat 
transfer from the gas medium across cylinder inner surface 
boundaries. Large variations in the maqnitude and direction 
of convective heat transfer exist from engine to engine and 
point to point within an individuial engine. The canditions 
influencing heat transfer are so complex and locally 
different through an engine cycle that rigorous mathematical 
interpretation is invalid. Radical Variations in 
instantaneous flow patterns allow gas velocities to be oanly 
Vaguely characterized. Non-uniformity of physical properties 
in the gas medium and surface materials exists. The 
combustion process Venice LCE sr iky difficult to model and 
further complicated by compression ignition. Two of the most 
serious obstacles question the fundamental assumption of 
fully developed, steady turbulent flow and the concept of a 
finite, instantaneous value of h. 

Tn operating engines, the entire cycle is completed 
within milliseconds, making accurate measurement of local 
instantaneous heat flux, gg, in itself quite difficult. 
Turbulent intensity in four stroke, CI engine cylinders is 


produced by the suction velocity of the induction stroke, the 


intensity being roughly proportional to it. To optimize 
fuel—-air micro-mixing, off-set induction ports may be 
incorporated to create swirling gas motion. In the 


compression stroke within a simple geometry cylinder, 





induction turbulence 15 suppressed. For swirl enhanced 
cylinders with combustion chambers cut into the piston crown, 
conservation of angular momentum in the cup created by 
viscous drag accelerates swirl to significant angular 
velocities as the piston reaches TDC. Radial "squish" 
velocities are created near TDC and combine with Swirl to 
produce a complicated toroidal motion. In the expansion 
stroke, rapid deacceleration of rotational velocities in 
direct injection cylinders are superimposed by ill-defined, 
small combustion induced velocities as hot combustion product 
gases expand; in indirect injection cylinders, prechambers 
deliberately channel and exaggerate the charge velocities to 
blow out of energy cells at levels up to S00QOm/sec C2]. In 
the exhaust stroke, the combuistion products are pushed out by 
the piston. The accumulation of these factors indicate that 
turbulent boundary layers are not steady, but are 
periodically created and destroyed in a manner beyond our 
ability to model. In addition, because of the thermal 
capacity of the gases, there is a thermodynamic phase lag 
markedly apparent in reciprocating engines; AT and gq do not 
pass through zero at the same instant indicating values of 
bh = to and highly fluctuating values immediately about it. 
Recognizing these obstacles, investigators have 
concluded that within obtainable accuracy it is still Best to 
model heat transfer aS a quasi-steady process and to 
correlate cyclic observations in terms of conventional Nu, 


Re, Fr 6ntumbers. Values of the various constants and 





exponents are determined in individual engine tests; a wealth 
of experimental data provides quidance for constant (s) 
selection for operating engines or design studies of similar 
construction. Hohenberg states [4, pg. 278391]: 
"Tt would be erroneous to look for a formulation 
which includes all possible influences on the 
heat transfer. This would only render the 
calculation more complicated rather than more 
accurate. The aim...must be to find out and 
mathematically describe the really essential 
factors influencing heat transfer in the engine." 
If the time for change in bulk stream conditions to diffuse 
through a boundary layer is small relative to other 
Significant times in the engine cycle, then boundary layers 
may be considered quasi-steady and at any instant be that 
associated with existing conditions. le Feuvre et al (C5) 


Rs 
ade 


approximating diffusing time as 6"/¥v,. and modeling gas motion 
in solid body rotation, successfully performed such 
calculations in a direct injection, swirl enhanced, CI engine 
Operating at high RPM (2000 RFM. As for the existence of 
phase lag, Annand (6] states that the actual quantity of heat 
transferred during this period is minute; the error arising 
from Bee use of a finite h at each instant of the piston 
stroke 1s negligible. 

The majority of the heat transferred within Aa 
reciprocating ak engine is due to turbulent forced 
convection. On the basis of experimental data and theoretical 
considerations, this - relationship 18s expressed by the 


formulas 


Nu = constant (R,)" Fy (14) 


i 





Recognizing that Fr = 0.7 for common gas mediums due to 
predominance of Q., and N-, substitution of the ideal gas law 


identity, p = F/RIg, into an expanded equation (16) yields: 
foes COmstant ere nar ee C7) 


Assuming that actual gas velocities will be at least 
approximately proportional to mean piston speed, ae function 
of Sp is used with the multiplicity constant determined by 
engine design. Tg is used vice the average, (Tg + Tw)/2, as 

the calculation of Tg from measured instantanecus pressure 
and gas law application already represents an average tor the 
cylinder contents. A consideration of gas disassociation 15s 
neglected, as any attempt to quantify it would be eclipsed by 
the previous assumption af uniform properties. During the 
combustion process, gas product composition is assumed to 
equal the burned composition. Tables and Figures of thermal 


properties (uy, C Fr, §) for individual gases and CHa, lean 


p? 
hydrocarbon fuel combustion products are provided in Appendix 
C 6,71]. Cyclic local instantaneous surface temperature 
profiles, from which actual measurement of q 1S generally 
derived, are Remon aated by Oguri et al C8] toa show 
transient magnitudes typically less tnan SO°F for sooty 
surfaces and less than 20°F for clean surfaces. With the 
exception of exhaust valve/port entries which may typically 
oe temperatures greater than 1000°F, the temperature range 


exhibited by cylinder sleeve, head, and piston crowns are 


fairly consistent allowing assumption of a mean Tw. 


—— =» 
. ° 





Se WOSCHNI‘'S CORRELATION FOR MEAN INSTANTANEOUS HEAT 
TRANSFER COEFFICIENT, h_,:-IN CI ENGINES 
A recognized industry standard for four stroke, direct 


injection, quiescent, CI engine heat transfer calculation is 


Woschni’s correlation, stated in U.S. units as follows C91]: 


: 


a. : on = = ’ _- ~S,-.27 7-2. 95 
ie 15.480 F CC, SptCa(VIg /FV_) CF Fo) J" ~B ike 5 bts) 


where AL = ETU/hr-€t<—-°R 


es vas Tg. = instantaneous values at ignition start 
Fo = motoring pressure 
Based on suppressed swirl engine designs, the values of 


constants C, and C. are: 


1. induction stroke C, = 6.18 Cy, = O (18a) 
2. compression stroke C, = £.2£8 Cyr = 9 (186) 
Se expansion stroke (motoring) Cy = 2.28 
(18c) 
Cy = 90 
4. combustion-expansion stroke Cy = 2.293 
ail (18d) 
C. = 5.7054 x 10 “ £t/sec-°R 
Ue exhaust stroke Cy = 6.18 C- = Oo (18e) 
From equation (17), an exponent of 0.8 was determined by 


Woschni by successively changing one parameter at a time over 
maximum possible range. A consensus opinion of all recent 
published correlations is that 9.7 =n = 0.8, rertlecting 
agreement with the laws of similarity. Woschni developed the 
motoring/firing correlation and scavenging correlation on 


suppressed swirl design cylinders stating [(9,p.307114: 





"The gas motion generated during the induction 
phase dies away with time, owing to internal 
Pie bet i OM. In the compression stroke, however, 
dependent upon the piston or cylinder cover 
shape, an additional gas motion can  6b6e 
induced Csquish-swirld] whose intensity is also 
proportional te) piston speed. In this 
particular case, the engine involved had 
direct injection and straight forward 
combustion chamber walls, so that decay due to 
internal friction prevailed." 


A ig °° exponent is due to Woschni’s assumption of air and 
Gas product sa(Tq) « Tg*°and (TQ) « Tg’? at n exponent = 0.8 
(see equation (17)). TO approximate ill-defined velocity 


increase due to combustion, the difference F-Fo which, unlike 
To-To, sees a rapid decay upon combustion completion, is 
used, multiplied by instantaneous cylinder volume, V, and 


related to the cnarge mass at Beginning of combustion, 


a mew / RIG. - 
It is at this point that Woschni makes aA critical 
Simplifying assumption. Recognizing the siqnificant effect 


Of luminous flame radiation heat transfer in CI engines, 


Woschni arbitrarily accounts for combustion radiant heat 


transfer by affixing the value C. to equate actual measured C) 
integrated over the stroke with @ combustion determined by 


the correlation. Woschni states (€9,p.20735]: 


"There is no doubt that such an influence occurs 
in the engine, and it ais taken into 
consideration in lumped form in the combustion 
term as defined here, but it seems impossible to 
the author to separate it from the convective 
term. The only question that remains to be 
considered, therefore, is the error made regard 
to the tendency...for the radiation." 





In essence, for application to the New Cycle Engine 
expansion cylinder operating at significantly lower peak 
pressure conditions, laws Of Similarity allow significant 
extrapolation of Woschni’s correlation for forced convection 
heat transfer only. The correlation must be corrected to 
\ 


account + or the fact that radiant heat transfer is not a 


function of (F-Fo)°®. 


4. CALCULATION OF ORDER OF MAGNITUDE REDUCTION IN 


COMBUSTION FORCED CONVECTION HEAT FLUX 


The representative firing cycle F-Tgq versus crank angle 
plot published by Woschni for natural aspirating operation is 
shown as Figure 2-1 C9]. The range of additional principle 
operating parameters over which the motoring/tiring 
correlations were developed well representive the essential 
parameters of New Cycle Engine (NCE) expansion cylinder 


operation [9]: 


Faramneter Woshni (W) engine NCE Expansion Cylinder 


REM 900 max Ser 
Sp 9m/s max C29.6 f/s] 33.00 ft/sec 
supercharge = Max 1.0 
pressure 
ratio 
Tw nominal So00°K [(900°R] FAO R 
ie, eS unknown 1.22204 Ft 


Bore and stroke dimensions were not provided, however from 
known maximum data point values Sp = 9 m/sec and Neey = 700, 


the approximate stroke length L may be determined: 





= 4&0 Some eGo) (So) eas m LO.98 tI 
Ne cong (22) (900) approx. 


For conventional CI engine geometry, bore diameter B is 


approximately equal to L; therefore for effective comparison 


of Hh, proportional to B°“ in forced convection heat 
transfer, consideration Of bore diameters of Similar 
magnitude dimensions may be neglected. The significant 


difference between combustion-expansion stroke protiles are: 
1. An order of 15 magnitude increase in peak 
combustion pressure in Woschni’s 


conventional high r CI engine. 


c? 
EF A reduced Tg decay rate, thus increased 
(Tg-Tw) forced convective heat transfer 
potential in the New Cycle Engine 


expansion cylinder. 


Sie A delayed fuel injection/ignition start 


ia the New Cycle Engine xpansion 

cylinder. 
Assuming that the magnitude of mean instantaneous gas 
velocities through the expansion strokes are roughly equal, 
the approximate ratio of forced convective heat transfer loss 
in the respective engines may by obtained by determining the 
mean “ratio of instantaneous heat flux, gq, over a YO? crank 
angle domain from ignition start. A selection of 90" 1s 


based on: 


ee Establishing a normalized basis in which to 


compare the order of magnitude reduction in 





forced convection versus radiant heat 
transter 5 the crank angle domain over which 
luminous flame radiation 15 significant in 
medium speed CI engines less than 1000 REM 
1s nominally 70° from combustion start C8Sil. 

ne At A the slope of pressure and 
temperature decay transients induced by 
combustion in both engines become small. 

Setting the radiant heat transfer constant Cy = OO in egquatian 


(18d): 


Gu/Fnce & Ag (Ta-Tw) y/A, (Tg-Tw) ine 


CFS 8Sp°8 ye “Tg* TIL Tg-Twiy 
eh ia co = << 
CPS 8Sp° 8 yet “Tgt FIC Tg-TwI yee 


Neglecting & consideration and assuming equal magnitude gas 


velocities: 


oh cP ®8/Tg*71LTg-Twl, 
Cx a= sso. 2. Ore 
+S a Ss wed at == ep 
Bee REO? / ign? Ct 1G = TI yan 





where Tw = 900°R 


Table 2-1 summarizes the calculations based on Figure 2-i and 
Appendix I data at 3S? increments from combustion start. The 
mean heat flux reduction ratio may then be determined by 


numerical integration using Simpson’s rules 








= _ 9O : ae 
INCE Forced convection : Q 

= ee qw qw qw qw Tu 

— Cl—— +4 +2-—— +... ¢2—— $4 +1— 
: : Snce  INCeE_ INCE. . Wnce_ 4nce__INCE_. 
Ow (yo =? 10° B80) So 909 
ANCE ATO" 

= AS°®/3 0328.67] = 6.987 approx. 
Soe 





Pectul ation 


forced 
Engine 


increments 


Crank angle 


past ignition 


start 


20 


29 


50 


39 


40 


45 


00 


59 


60 


69 


70 


79 


80 


85 


90 


p 


ko/cae 


40 
60 
63 
40 
55 
47 
40 
33 
26 
22.5 
20 
17 
15 
12.5 
1 
10 

; 


8.9 


order of magnitude reduction 
convection 
expans1lon 


Wosfini engine (W) 


Tq 
1K 
875 
1400 
1525 
1600 
1650 
1430 
1600 
1540 
1480 
1440 
1400 
1365 
1330 
1305 
1280 
1255 
1230 
1215 


1200 


cylinder 


P 


ata 


33.71 


08.07 


60.97 


58.07 


JS. 24 


45.49 


38.7! 


31,94 


25.16 


21.18 


19.36 


16.45 


14.92 


TABLE 


Ty 
OR 
1575.0 
2520, 0 
2745.0 
2880.0 
2970.0 
2934.0 
2980.0 
2772.0 
2664.0 
2592.0 
2520.0 
2430.0 
2394.0 
2349.0 
2304.0 
2259.0 
2214.0 
2187.0 


2160.0 


NCE Expansion Cylinder 


P 


ata 


4.49 


4.20 


4.09 


1,98 


1.49 


1.42 


1.33 


1.3 


1,29 


= 0 


at 


iG 


Tq 
7 
1896.6 
2252.6 
2571.2 
2864.1 
3134.0 
3064.2 
2998.7 
2939.1 
2884.8 
2835.2 
2790.0 
2748.8 
2711.3 
2677.1 
2646.0 
2617.8 
2592.3 
2569.2 


2548.4 


Crank 


combustion 
q from Woshni’s test engine to New 
operation 


Cycle 
angle 


Wy! SNcE 
equation (19) 
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FIGURE 2-1: Woschni's test engine firing cycles natural 
aspiration L9 
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eo . ESTIMATION OF MEAN GAS VELOCITY MAGNITUDE ,U 
ei. CONSTRUCTION OF CHARGE IMPULSE VELOCITY ,Ui, TRANSIENT 
Acceptance of the general applicability ofr Woschni’‘s 
correlation to the New Cycle Engine expansion cylinder 
balances on the assumption of comparable gqas velocity 
magnitudes. While geometric flow patterns are indeterminate 
in preliminary design, the magnitude of the charge impulse 
into the expansion cylinder leaving the regenerator port can 
be ascertained. Recognizing the inevitable deacceleration of 
this aimpulse in a simple geometry cylinder on the expansion 
stroke, aA conservative estimation of the mean Bulk gas 
Velocity at the walls, Vexp, through the expansion stroke can 
be made with minimum error. Substituting Vexp directly into 
the steady, fully developed turbulent flow equations (14) and 
(15) for flat plate and pipe flow should, if the magnitudes 
of cylinder velocities are in fact Similar, be reasonably 
Glose to h, calculated from Woschni’s correlation. To reduce 
uncertainties, the relationships will be first compared in 
fieaeGgeing cCONditions from TDC to 45°, prior to ignition start. 
Dent et al CLIO] successfully used a similar method in 
approximating forced convective heat transfer loss ina four 
stroke, direct injection, high swirl, CI engine operating at 
peak TDC Swirl ratios equal to 20, Fredictions using 
Woschni’s motoring correlations failed uniformly low compared 
to measured, spatially averaged, instantaneous heat flux. 
This 15 expected, considering the central assumption ef 


suppressed swirl inherent in application of Woschni’s 





correlation. By substituting directly into the flat plate h 
equation (14), measured instantaneous angular gas velocity w, 


at i equal to bore circumference (2rr), remarkable agreement 
was obtained. Agreement was not as good at firing conditions 
due to the transient combustion processes. 


From ideal gas law, at a given angle 9 in the expansion 


stroke, charge mass, m, in the expansion cylinder: 


EV (2116.2724 lbf/tt“—-atm) FY 
Lemiim: Ts ee 
RTG (55.90) Faq 
= 39,262 FV (24) 
g 


Taking S° crank angle as the unit increment of time at 818.3 
REM; the charge mass rate, entering the expansion cylinder 


1s: 


' Am (5° ) 
AG Gae.) (Neem / OO) 260° 
= fm (5°? ) (22 
~9010184 sec 
The continuity equation for impulse flow leaving the 
regenerator port cross-section area, A,, equal to 0.80 £t= 


1s: 


a) 


Therefore, By rearranging equation (23), impulse velocity 


leaving the regenerator port, vi, can be determined: 





maRTg Me, (53.90) Tg 


Vietysec a a =. 5 oe 
F'A (2116.224 Lbt/tt7-atme¢.080) 


2 21857 m., Tg 
Se aa a (24) 
EF 
A summary of the calculations using instantaneous F, Tg, V 
Values from Appendix A at S® increments is provided in Table 


°- 
we 


Be Charge mass as a function of expansion stroke crank 
angle is plotted in Figure 2-2, indicating that the main 
impulse from the compression cylinder via the regenerator 


x1it port completes 10° past ignition start. 


At 45° combustion start “%mo..) = 2. 1905x104 = 89.5% 
2. 4527x10 - 


At 3o° main impulse completion - 
AM ot a] = Ho 28874 10" Ss as /a 
ee F4557K10 ™ 
Ui aS a function of expansion stroke crank angle is’) plotted 


im Figure 2-3, indicating a peak Ui equal to 3519.95 ft/sec 


mear fuel injection start (20°). 


44, 





a es Sg 


+ 
Ci 


Calculation of charge Ui leaving the regenerator port into 
the New Cycle Engine expansion cylinder at S° crank angle 
increments 
Crank Angle m Ameo m Ui 
q° lbm ATane—= lbm/sec ft/sec 
eq. (21) eq. (22) eq. (24) 
TDC o o o . 
o 1. 5550E-4 OAL oo (Nap 7 41.82 
1 7a /OFOE-4 -6154 26043 132.14 
a 1, PADS Be 1.260 ae Gs aes 
20 4,50821E-3 2. 45 oe 4069 Je, 
ae Se ke any aot) aoe dics 438.13 
Aw 1p oe 4.625 4.5413 Sioa og 
sad eee 4.437 4.357 AG a Aa 
40 anes Le ae 220614 ea 
435 eee 1.414 oo (ize, 
wt) Ageoo7 e— 2 G29 S816 «7641 degra fo 
oo Soe a QO. 9028 SOO 7S ees 
50 Zee oe 0.90426 ees Fe 
65 2. 29923E-2 0.0596 0.05852 ees 
70 a noeosoe a ~9641 0.96294 eee ys 
To aay Ei Ome2a lt) 0. 2465 81.44 
80 oa 2 0.2049 Ono) 2 2.76 
85 Sooo On = O,.1687 Toa / Sa 
PO eo eoe 2 Ones Oo ye Obs 
oo ees 7 | O.1167 ~ 1146 yt (icy 
LOO ee) Oe worst Ore © acs 46.356 
10S ome oe 2 a OB 7 mma l ai 41.354 
110 <.4188E-2 0.0702 O,06393 owe bs 
115 on 4247E-2 Ona 7 © Ooo ets saa) 
120 fe eo Se O,OS 11 O,05018 =9.84 
| oad ae 4i2423E-2 0. O436 O,O04281 eee ic 
LO 2a? gas ~I9372 en sw OF 
a 2.4411E-2 ro O.O8 142 =O. 41 
140 2.244 28E—-2 290270 OeiZzoo lt Laon 
1435 2.24462E-2 ass 2. O2247 ea we. 7 
1350 =. 4480E-2 I. 0185 POLS 7 12.48 
So =.4496E-2 V.OLS3 OLSo2 10.94 
160 Sa toOF ea ~O1LE0 O.O1L2E7S Fake 
Loo eat Ge O.Q0101 0, 00992 Tee 
170 apenas 7 te a O.GO74 O20 727 Saas 
1735 Zatonae— 2 0.0031 0.90301 Seog 
BDC ooo oe 0.0032 QO,Q0311 a 7 O 
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IGUR# ¢-2: Expansion cylinder mass charge, m, versus 
expansion stroke crank angle 
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ie < VERIFICATION OF TURBULENT FLOW 


To 


verify turbulent flow conditions leaving the 


regenerator port, calcttllated Re at 20 using Te equal to the 


regenerator 


Me7o inch iss: 


R@=¢)0 











port cross-section minor dimension (depth) of 
Uil- Vil *Uil-.. i 
= = — eee (25) 
vy M/p M RTg 


From Appendix C, Table C-1, w(Tq) and C,(T) for air ist 


wm lbm/ft-sec = 2.163x%1077 (Tg) *O** (26) 
(C1) 
C, BTU/ibm-hr = 0.224 + Pee oetoman ia Cup 
dee) 


Substituting 


instantaneous values of F, Tq from Appendix A 


Seemed) identity equation (26) into equation (25): 


wn) 9 


RE=(0 


6.10171 atm 


2194. 29°R 


CS oe. 7 ay Ley Cet i6. 24 lbf/#t*-atm) (6.10171) 


(2.162x1077) (2194, 29) O49 S%.90) (2194.29) 


a 
1.1476x10 turbulent 
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ec ESTIMATION OF EXFANSION STROKE ,Uexp 

The Calculation of the mean impulse velocity, Ui, for 
Biempuroose Of Vexp approximation, is correctly restricted to 
the main impulse domain TDC to SS? only, neglecting residual 
1OA charge mass regenerator bleed. Combustion induced 
perturbations during crank angles 42°-S335° may be neglected. 


Ui may then be calculated using numerical integration by 





trapezoidal rule and 1s equal to: 


Ui ay = gal [99° Wi de (28) 
SoceASS | a 


(AS°/2)0 Vigne tise +2Vizgo+-++WVisgetisge] 


ASS ° 


(AS°/2) 05979.03] = 271.77 ft/sec 

oso? 

A reasonable assumption must be made as to the expected 
value of UVexp over the entire stroke. The principle factors 
causing charge velocity deacceleration include cross-section 
area expansion exiting the regenerator port, internal 
friction, cylinder volumetric expansion, and angular momentum 
deacceleration of any existing rotating flow patterns. 
References (5, 10, 11, 12, 13] unanimously agree that in swirl 
enhanced engines, peak instantaneous swirl at TDC 15 much 
Greater than the mean swirl valves over the stroke. 


fe Feuvre et al £€3S] calculated for a 4.35 inch bore direct 


injection, CI engine operating at 2000 RFM, that 104 velocity 





decrease occurred each stroke due to viscous dissipation 
alone. A conservative estimate of Vexp = (Uitnpe-sso)/2 is 
therefore assumed; selection of this value should reduce any 
Peo in either direction. nis Horeca that Vitpc-rpe fram 
Table 2-2 im firing conditions equals 1092.9 ft/sec. 
Combustion induced velocities due to charge gas expansion in 
Poemeeeetagmition (SI) engines have been calculated by Amnand 
C2J to be on the order of 32 ft/sec, with expected velocities 
im CI engines less due to diffusion flame configuration. 
Thus, let Uexp = (Uitperisso)/2 lee Cts Coe 
6. COMFARISON OF he CALCULATED FROM WOSHNI’'S CORRELATION 
AND THEORETICAL MODELS 
5.1 EXPANSION STROREs; MOTORING CONDITION 
Using the above first approximation of Uexp, estimates 
of ae from theoretical equations (14) and (15) for steady, 
fully developed flow may be reached. Expanding pipe flow 


equation (15) at 1. = B and U = Uexp yields: 


Ap = -O23(K/B) Rep Ser! 44 = .0235(K/B) CUexpB/y,2° FPri/9 


Substituting the identities vy, = p/p, K = CHAP rE and 


p = P/RTg: 


ed 


a 


p = -O23C, actPUexp/aRTg1° (Oyo ep, ee. 


Substituting the C567) Hidde eedentitlies from equations 


(27) and (28) and appropriate conversions: 


ei) 





hp = -O230. 224+. 262x107 7 Tg Il 2. 163x107 7 (Tg) * O44 (3600sec/hr) 1# 


C(2116.224 1lbf/ft*-atm F Uexp)/(2.165x1077Tg* oF 9R TQ) 1° Sx 


a7 =p, 72/5 


which simplifies at known R, Vexp, Fr values to: 


CSG) at) is 
BS, ME ac eat Geax Om ag? —___. (29) 
ETU/hr-#t*—°R Be! Tg: o7! 


An identical expansion of the flat plate equation (14) at 


ee = wH, expansion cylinder circumference, yields: 
(135.9 p)°1S 


Tg] ——______ (230) 
Re Tg o7t 


My _ = 11.674023.24+2.62x% 1077 
Pe 7 hr —tt~—° RK Ee 

A summary of A calculations for the charge induction 
monitoring angles, TDC-43°, comparing equations (29) and (320) 
to Woshni’s motoring correlation equation (18c), is provided 
arable 2-3. The agreement 1S quite qood, even with 


acceptance of reasonable error in VUexp approximation. 





Lane 


if 


Comparison of h, calculated by Woshni’s correlation = and 


theoretical models during charge induction into the New Cycle 


Engine expansion cylinder 


Crank Angle Ne An hy 
1° eq. (18c) eq. (29) eq. (30) 
~ 16.994 tae eet 7s 
10 eae Lo 1 18.143 eee 1s 
cs 24.006 le Ss est fF aioe 1. 
20 meee oe la cF ae 496 
eo i la Seco 7 ro 7 Gna 
oO 24.226 ale) om 07 39. 4038 
ed oo. SBE 20.414 Sees 
40 Sloe a7 wos 7 ee 7 7 
435 27a Lod Be aes lee 7 


h in BTU/hr-ft<-°R 





a EXHAUST STROKE 

Attempts to directly estimate the mean bulk gas 
velocity at the walls during the exhaust stroke, Uexh, in 
light of preliminary design restrictions and post-combustion 
ambiguities, are not valid. 

Pmncatect method allowing h calculation is to correct 
equation (29) by a factor equivalent Co the 
scavenging/motoring Vexh correction observed by Woschni for 
suppressed swirl geometry, direct injection, CI engine 


cylinders. From equations 18a and i18c: 
Uexh correction factor = (6.18 Sp/2.28 Sp)°? = 2.22 


When the regenerator exit valve opens, the 
Significant flow pattern will be parallel to the bore axis, 
up the sleeve, and through the regenerator port; flow 
direction is therefore analogous to pipe flow. From Appendix 
C, the identities w(Tg), Cy (Tg) Pre (vy), and (Tg) for lean 


C Hae, hydrocarbon fuel combustion products at F = 1 atm, 


6 = 0.75 may be calculated: 


1.47x10°7(T7g)*% 1.47x1077(Tq)*7 


& Llbm/ft—-sec a 
1 + .0274 Dee) 7 (57) ) 


= 1.441x%10°7(Tg)*7 (31) 

(c=) 

C_ BTU/lbm-°R = 0.247 + 0.243x1074 Tg (32) 
P = 
Cea) 

mar = 05 + 4.2(7-1) — 6.7 (3-1) Ges) 
(C4) 


5 (Tq) 


fel 


Appendix C, Figure C-1l (34) 
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Applying the Uexh correction factor = 2.22 and the identities 
pectg) , C4 «Tg? defined above to pipe flow equation (29) 


yields: 


hp = 2-22(.023) (247+. 243x107 "Tg) (1.441410 7Tg* 7) (3600sec/hr) # 


C(2116.224 lbf/ft*-atm PF Uexp)/(1.441x1077Tq* “RTg) 3° Fx 


which simplifies at known R_); Uexp values to: 
1 al Sd eee 
An - MGs soCendgeer 2. 45x 10n1g] —— pp 2/s 
BTU/hr—-ft*-°R B°“Tqg* 4 
(55) 


From Appendix A; F = 1.95 atm and 2344.59°R = Tq ee a ec Ex 
Loe EDC 


Baeotam the exhaust stroke. & comparison of A calculated 
from equation (35) and Woschni’s exhaust stroke correlation 


equation (186) follows: 


crank angle Tq SCG? Fr cy) he An 
OR (°K) eq. (34) eq. (35) eq. (18e) eq. (25) 
RDC 2439.06 Las o 0.698 17s oo7 17 2Oo£e 
Cle oa) 
TDC 2244.58 1, 230 Cet Pat 7 17 ok 


ae) 
A in BTU/hr-ft*-°R 
The above series of calculations suggest that to a first 
approximation, application of Woschni’s correlation in 
determining forced convective heat transfer loss for the New 


Cycle Engine expansion is valid. 





ALCULATION OF UNCORRECTED HEAT TRANSFER 


LOSS/HEAT RELEASE FER CYCLE, ft 


Table eo summarizes the calculation Of total 
instantaneous Neat transter rate, a, through a complete 240° 
firing cycle of the Néw Cycle Engine expansion cylinder using 
Woschni’s correlation for h.. For calculation of AL in 


combustion, motoring ae is generated graphically in Appendix 


A, Figure A-1l: 
CeTU/ssec = n Any CT g-tw) Ct he 73600 sec) (36) 


Due to the establishment in preliminary design of a simple 
cylinder geometry, Aut as a function of crank angle may be 
determined directly from instantaneous V values; defining 
instantaneous gas volume height as L, and neglecting cylinder 


sleeve clearance: 


V = WE L, 
4 


Au, = 20E“ + wBL, = 2B“ + wDC4V/qB*] 


sh 4 
which at B= 1.22204 ft: 


Aut = 2.246 + 3.27S5CV] Cav7) 


- 


#t~ 


The distribution of ie and @Q as a function of crank angle is 
plotted in Figure 3-1. It 1s immediately clear studying the 
Meecastribution aver the combustion stroke SS° to BDC, that 


= 


Woschni’s correlation fails to correctly extrapolate luminous 


Ci 
in 





flame radiant heat transfer loss to low Fr operating 


c 
Samedi tions. No increase in h, occurs due to the predominant 
influence of increased temperature: associated qas layer p 
decrease and uw increase overcome the effect of E increase, 


thus lowering 5h modeled by laws of similarity solely an 


ea? 
forced convective heat transfer. The low r_ and delayed fuel 
injection operating coanditions effectively eliminate any 
remarkable pressure rise in the cylinder, already at 45° 
down-stroke at ignition start. Recalling Woschni’s coupling 
Of radiant heat transfer in combustion to C5 (PF), the 
desired “lumped form" term eftectively disappears. 
Uncorrected total heat transfer loss per cycle, le, is 
determined oy numerical integration under the O curve using 
Simpson's rule: 
ae — 00101840 Bppct Reed 20h eee, eas - 


BTU/cycle = ps 
FeO FQ 0 t4Oy geo the ned (28) 


= .901018407141.612)] = 2.42454 BTU/cycle 


% 
% 


From Table 2-2, total charge mass per cycle equals .024527 
lBbm. At known specific heat release per cycle, A heat, equal 
to SOQ BTU/1lbm charge, the ratio of total heat transfer 
loss/heat release per cycle, hy, equals: 


R+ = =.242424 BTU 
uncorrected (SOQQRTU/I1 Dm) (.9024527 lbm/cycle) 


= 19.7617 


ey 
{}- 





The inability of Woshni’s correlation to account for radiant 
heat transfer in low r_ operation is Clearly demonstrated by 
repeating the numerical integration with the combustion term 
in equation (18d) eliminated (C., = ©): 

Simpson's Multiplier = 7038.088 

QL = 2.28972 BTU/cycle 


Ry = 19.474% 
uncorrected 


Considering the accuracy of the correlation in general, for 
preliminary design prediction purposes setting Ce = © 415 
considered valid. 

To correct for luminous flame radiation heat transter 
loss, a value of qa7aeee radiation must be determined which, 
in comparison of auvance forced convection calculated in 
Chapter Two, will allow first approximation oft the required 
Ry correction factor. As there exists no theoretical basis 
for couplinq magnitude of luminous flame radiatican to 
(PoP) 8, an attempt to determine the desired ratio from 
(P-Po) profiles must be disregarded in favor of calculation 


founded on known principles of radiation heat transfer in 


diffusion flame combustion. 


Cf 
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Bateulation of 


» 


cycle @ using 


increments 


crank angle Re Aut 
a BTU ft* 
hr-#t7—°R 
eq(18) eg (a7) 
charge induction equation (18c) 
is 16.944 eee 7 
10 oO ako Ze) 
15 24.906 2. 445 
20) ee Oo Beenie 
eae oO, aes 
30 4. 226 ee 724 
oo 2. 888 Z2eGe7 
40 51.168 Sorenery 
45 27 AauGs eee t 
combustion—-expansion eguation (18d) 
ow 26.086 Joao 
Sle ee OF oe I4B 
60 22.966 52 741 
oa 20. oa =e P44 
70 Logan @ 4.1352 


TABLE 2-1 


New Cycle Engine expansion 


Woschni’s correlation 


cn 
G 


Fn Red 
Ta-Tw 
OR 


Tw=900°R 


1249.04 


Deeg el 


1295.75 


1329. 26 


beg. 7 I 


1294.29 


2 ow ot 


ee orcs 7. 
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BTU 


Sec 


eq (246) 


12.856 
Nowe? S 
ole LES 
202190 

Deere 
oe Oe 
ae BI4 
a7 Oo 


seals 924 


ae OBB 
=F. 445 
435.028 
ye sc 


46.700 


firing 


angle 


crank 


angle 


oO 
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100 
noe 


110 


hese SA 


Ne 


RTU 


hr-#t*-9R 


17.084 


is. /42 


14.632 


13.713 


iz.841 


Lege 


11.457 


10.8635 


10.428 


Loni 


8.830 


8.632 


8.467 


8. 3.34 


8.229 


6.1951 
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8.044 
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(CONT ° 


ae 
ie 


5.927 


S405 7 


&./11 
Cee? 7 
6.8/1 
6. F730 
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7 | 
72051 


Oe) 


Tg—-Tw 
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2059, 11 
1984.75 
1 38 
1889.78 
1848.18 
Po tices 
7 oy 
1746.03 
1717.83 
1692.29 
WeG7s22 
1648.44 
1627.32 


Lol eee | 


BTU 
Sec 


45.444 
40.775 
38.584 
34.789 
25.002 
23.478 
32.148 
20.877 
at), OOSG 
29.089 


Sonal 


24.917 
24, 26.5 


24.260 





TABBE S-1 
crank angle is 
a BTU 
hr-ft=-eR 
exhaust equation (18e) 
EDC 17 .oev 
185 17.604 
190 172 OL7 
ey 17.624 
200 17.4549 
20S 17.664 
210 17 wos 
eo eo 
2eQ) 177 OF 
a 17.7 
ao 1/7 a7 So 
Zoe 17.748 
240 7a ol 
245 17.774 
es) 7 So 
oS5 Iya oa 
260 Lo. 
Zoo eee el 
70 172 oo 
oa) 17.842 
=80 La ee 
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ee 
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pee ok 
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LSGye 7 7 


1504.41 


1500.90 


1497.47 


1497.47 
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BTU 
sec 
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49.035 
48.121 
Aes 
46.037 
44.870 
43.628 
42.3508 
40.929 
39.496 
38.012 


26.408 
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FIGURE 3-1: Expansion cylinder firing cycle he and Q versus 
crank angle; Woschni's coorelation 





CHAF TER FOUR 


RADIANT HEAT TRANSFER IN THE NEW CYCLE ENGINE 


EXPANSION CYLINDER 


ie THEORETICAL BACKGROUND 
When a temperature difference exists between two 
bodies, part of the thermal energy is transferred as 


electromagnetic waves of radiation at varying wavelengths, A, 
characteristic of the body temperature and identity. In the 
electromagnetic spectrum, thermal energy 15 propagated 
through .OL wm = A £ 100 pm. Bodies above absolute zero 
temperature continuously emit energy which, when striking 
another body, may be partly absorbed, &%, reflected, %§, or 


transmitted , 3, where: 
Bes sce ee = (55) 


Thermal radiation reflected from and transmitted through a 
body 1S again absorbed, reflected by, or transmitted through 
ray additional bodies ultimately establishing an equilibrium 
state. An ideal body that absorbs all radiation is known as 
a Black body; if no radiation is transmitted, as in the case 
for most industrial liquids and materials, the Body is known 
as Opaque. Combustion chambers covered with a thin film of 
Si] contaminated with soot particles will absorb A 
Significant percentage of incident radiation, FT = Q, 

Thermal energy emitted from a body unit surtace per 


unit time at a particular wavelength, the monochromatic 





emlisslve power, Eas can be determined aS aA FUNCTION oF 


temperature TCP RF) by Flanck’s Law: 


E. _ = 2mee,*\ "Cee, /eSb*!-137! (40) 
erg/sec-cm~ 


speed of light 2.997902x10!9 cm/sec 


where: e = 
Cy = Flanck’s constant Be 7 erg-sec 
Ch = Boltzmann's constant 1.38026x10 °° erg/? kK 
A = wavelength in cm 
T = temperature in °F 


ES will auncrease with increasing A to a maximum value then 
decreases, the distribution shifting toward the shorter 
wavelengths at higher body temperatures. The relationship 
between 1T(°R) and AL, is determined by Wein’s Law of 


Displacement: 


Ea T(PR)D = 2898 cm-°F 
The total radiant emissive power, E—, may therefore be 


calculated at a given body temperature by integrating the 


area under the curve: 


E= (["c&. aa (41) 


For a given Black body, substitution of equation (40) into 


equation (41) yields tne Stephan Boltzmann Law: 


Z 4 7 
Eplack Body = %o! (oR) (42) 


BTU/hr —-ft* 


where: o«, = Stephan Boltzmann constant 1.713%107 7BTU/hr-Ft=°R* 


&4 





Real bodies, however, are not absolutely black. Those 
bodies of a continuous radiating electromagnetic spectrum 
where the E. curve 1S proportional to that of an absolute 
black Body at each A are known as grey bodies. The ratio of 
emissive power of a grey body to an absolute black body is 
known as the coefficient of emissivity, €: 


eagrey body (7) = € (4.2) 
E black body (T) 


From an analysis of equilibrium emittance and absorptance 
between a grey body and its Surroundings at the trivial case 


of each equal to the same temperature, it can be demonstrated 


that © ais equal to jo (Kirchoff’s Law). This important 
phenomenon, strictly correct only within the defined 
limitations, is fortunately observed to be approximately 


correct in many real conditions where temperature differences 
are not severe. Additionally, many real industrial materials 
radiate near grey at very high values of W&; instantaneous 
radiant heat transfer rate, Oe, 5 between two near-grey bodies 
land 2 of different temperatures at a distance apart may 


therefore be calculated from equation (44): 


: _ A ae 
R1,2 = Fy ar K 1 287 (Ty ~Ta 9) (2 br 73600 sec) 
Rtu/sec 
(44) 
1 
eee: £, 5 = Lri- i 
j= a — 
1 §2 
Fie = geometric view factor between bodies i and 2 





Tables providing formulas for € and FF calculation 
between multiple bodies at various common geometries are 
published in standard heat transfer texts. For gas oar 
incandescent solid particles radiating as a defined unit body 
within an encompassing simple geometry combustion cylinder, F 
may be pee tned equal to unity. 
es NON-LUMINOUS GAS RADIATION 
Zl CHARACTERISTIC HEAT TRANSFER EQUATIONS 

Unlike solid bodies which radiate across continuous 
wavelength spectrums, qases emit and absorb energy at direct 
wavelengths only in line spectra. For simple monatomic and 
diatomic gases such as Oo,No5 and H,, the line spectra are 
negligible and may be considered transparent (7 = 1). Por 
complex gases such as generated in exhaust products, however, 
radiant capacity may be significant. CO. and HO vapor 
predominate as radiant emitting gases in combustion products; 


the significant bands of radiation are C1i4]: 


CO, H.0 
Band 1 mes - 3$.02 pm eect - 3.27 AM 
Rand 2 4.91 - 4.80 ym 4.39 - 38.5 wm 
Rand 3 AS ai) ai ae Ly 12 oe yum 
In radiating gas, both emittance and absorptance exist 


throughout the gas volume. Fenetrating radiant energy will be 
reduced by absorption proportional to the number of molecules 
in the ray path as a direct function of individual gas 
partial pressure, temperature, and path length, Xx. This 


intensity attenuation is expressed by Reers Law: 


56 





= t-exp (tk, X) (435) 
where: Ey = monochromatic absorption coefficient per 
unit length 
A = radiant path length 
A gas volume is, Or CGOURSe., wmat grey. For ease of 


calculation an equivalent gas emissivity, €¢ 


g? is defined in 


accordance with equation (464): 


EG = i1- exp (Kile) (46) 


wheres k apparent gas absorption coefficient per unit 


g 
length 
le = effective mean beam length 
For a defined gas volume of V and boundary surface Airs Ll. 
may be determined from equation (47): 
Le = 2.44 (47) 


ft Aut 


® 


Ge between the gas and surrounding wall of emissivity GY MAY 
then be determined from equation (48): 


Q =€€ 


ie 4 (Sex es ae 
R ge wTofHTl Tg -Twti(t hr/S3600 sec) (48) 
BTU/sec 


For gas volumes containing more than one radiating gas, total 
Bg is calculated by addition of individual gas emissivities 


and an additional negative correction factor, —At - ete 


(a 
radiant band coincidence which reduces the sum. For internal 


combustion engine exhaust gas products from hydrocarbon 


6/ 





based fuels, the resultant total gas emissivity may be 


considered equal to: 


Ge =€ °+€& -a 
g co, HO 9g a 


Again, significant data base has allowed development of 
graphs for € es 


EG. HO’ 


Les and are published in standard heat transfer texts. 


and ~ OE as a function of FF, Tq and 


Empirical formulas based on experimental data for calculation 
Of emissive power of CO. and H»0 vapor useful for ai first 


approximation in U.S. units are C14]: 


pees avi@antO “(Pop 1,)'/*(Tg/io0)=> (50) 


Ws 


ETU/sec-ft= 


van 


4 Sat teva g = Sai : @ S a & = “40 = e C) 
Ey o a 1 D 67.28% LO PH_o it FY (Tg / 1 QQ) 
< oe 


ETU/sec-¢t* 


tJ 
br) 


FREDICTION OF GAS RADIATION HEAT TRANSFER 

References (2, 6&6, 8, 7%, 10, 14, 15, 16] state that gas 
radiation 15 insignificant in calculating high pe CI engine 
heat transfer loss. To ascertain if separate consideration 
Of gas radiation in the low i New Cycle Engine expansion 
cylinder may be required, it 1s sufficient to estimate O. 
magnitude at maximum It. and Tg in the exhaust stroke and 
compare to Q calculated by Woschni’s correlation. For a 
typical diesel fuel of C/H mass ratio equal to 6 operating at 
lean $9, distribution of combustion product gas includes 


approximately LOZ CO. and 10% H.Q water vapor Ci4i). 
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Beapeulation of Oe. using equation (48) and (49) at Appendix A 


parameter values at exhaust stroke BDC tollows: 


Q Woschni = S2.909 BTU/sec 


V = 1,43333 £t7 

Aut = 7.037 t= 

Tq = 2429.06°R 

P= 1.05 atm: Py g = 105 atm Pog = +105 atm 
ip = _3.4(1.43333) = .6925 ft 


72937 


Ag (1.7125¢1077) £60105) (26 BS) “edera@e, ~*~ C7masm 
= LO0 


= 3.604 ETU/sec 


(1.6738x 1077) (.105) °F ¢. 6925) °° (2439.06, 9° 9(7.037) 
a 100 


g 
ee 
O 

ee 


= 90.226 HTU/sec 


- 


Up. = &5.8:2:0 BTU/sec 


In the above calculations, corrections for DE as E and 


Ww? 


emissive power of the cylinder surface proportional to soTwe, 


all which would reduce the total, have been ignored. The 
estimated value equal to 32.83 BTU/sec is well absorbed into 
wee 709 EBTU/sec. Noting that for all other crank angles in 
the exhaust stroke both Tq an lp decrease, plus in the 
combustion-expansion stroke luminous flame radiation exists, 
ie may be assumed that separate consideration of gqas 
radiation in the New Cycle Engine expansion cylinder heat 


transfer calculation is not required. 
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=. LUMINOUS FLAME RADIATION 
te 1 FUNDAMENTAL CAUSE 

The phenomenon of flame radiation is quite complex, 
dependant intimately on the physical and chemical properties 
Of combustion, pressure, temperature, and a myriad of engine 
and fuel parameter variables. Two basic flame types are 
produced with respect to radiant capacity: non-luminous flame 
characteristic of nk engines, and luminous flane 
characteristic of CI engine diffusion flames. In SI engines, 
pre-mixed fuel and air are adiabatically heated in 
compression then iqnited by a spark resulting in combustion 
by turbulent propagation of flame through the charge from the 
point of ignition. The mixture Burned at the beginning of 
flame travel xpands, then recompresses, reaching 
temperatures Rhundreds of degrees higher than that which is 
burned last. This recompression of combustion product gases, 
particularly CO. , yields visible yellow radiation. No 
intermediate solid particles are formed, nOwever , and 
apparent emissivities of combustion product eee as a 
function of near transparent A bands is small. The driving 
potential for radiant heat transfer, proportional to lean 
diminishes more rapidly then forced convection heat transfer 
proportional to Ts the contribution of instantaneous radiant 
heat transfer, Qe, integrated over the engine cycle is 
insignificant. 

In CI engines, the fuel spray injecting into 


adiabatically compressed air near piston TDC is ignited by 


7O 





spontaneous explosion of a micro-mixture or fuel and 
entrained air, thereafter burning diffusively. There are 
many sources of ignition and flame tends to conform to the 
fuel spray shape until dispersed by gas motion. In the 
chemical combustion of hydrocarbon based fuels, cracking to 
simpler structures and polymerization occurs. When 
insufficient local 0. exists to feed the rapidly increasing 
reactivity of the mixture, thermal decomposition of C-H and 
C-C bonds occurs resulting in incandescent carbon soot 
formation radiating in a continuous spectrum. Alcock et al 
C11ljd, im a classic high speed photographics study of a Meurer 
direct injection, medium speed CI engine, clearly observes a 

luminous flame front filling the chamber within 20° = of 
Ignition start (10° BTDC), glowing white hot at peak flame 
temperatures equal to 4500°R, then cooling by expansion to 
dull red. Unlike SI engines where combustion completes early 
in the expansion stroke, the CI combustion process extends 
late into the @xpansion stroke. 

This difference between local QO, concentration = and 
cracking rate 1s the principle cause of soot formation = and 
luminous flame radiation. The magnitude of associated heat 
transfer canmot be ignored, and has been measured by Annand 
C3d, Sitkei C141], Dent et al C102, Flynn et al C151, and 
Oguri et al C8l as typically 20% to 30% of total heat 
transter loss, &) , In non-swirl enhanced four stroke, direct 
injection, high -r Cl engines under varying design and 


Cc 


operational conditions. Oguri et al (C84 demonstrate tnis in 


fi 





Figure 4-1, showing @./Q, Between two such different size Cl 
engines operating at varying RFM, load, and cylinder wall 
soot contamination: Engine Ais a S.97 inch bore natural 
aspirating engine; Engine Bis a le.2 inch bore supercharged 
engine. In very large diesel engines the contribution to Gl. 


Of radiant heat transfer 1s even qreaters J. Butler of 


Dexford Co., United Kingdom, estimated the radiant component 
of Ch. as 73% for a 40 inch bore diesel (81). In any case, 
measured peak flame emissivities, Oey approaches unity 


(1.9), and while reaction zone temperature is difficult to 
determine, a consensus opinion shared is that peak radiant 
temperature is much higher than Tg. Chapman et al state 
muy aeeoc Lsls 

"the apparent radiant temperature is much closer 

to the flame temperature than the average bulk 

gas temperattre." 
Feak flame temperature, Te, reached in high r_ CI engines is 


typically 42900 4 200°R. 


aeitgn KHAN ET AL COMBUSTION SOOT FRODUCTION MODEL 

While a valid quantitative calculation or soot 
concentration in New Cycle Engine combustion is not possible, 
a brief study of a soot production model for a simple design, 
direct injection, quiescent CI engine will illustrate the 
critical parameters on which soot formation is based and 
establish a foundation for empirical estimation of peak New 
Cycle Engine @,. ~ Khan et al [18] have developed 


semi-empirical equations for soot formation ina fuel jet 


mixing model applied successfully to a 9.93 cm bore engine at 


“J 
bd 
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20 40 60 . 80 
%, MAX, IMEP 


Q2/Q7p for two 4 stroke, direct injection CI 
engines under varying operating conditions as 
observed by Oguri et al 


2 


lOO 





2000 REM 


This simple model ignores complex stratification 


and consideration of fuel droplet evaporation; all unburnt 


fuel is assumed to be in the vapor state. 


injecting 


tae a cross-stream air tlow iss shown mn 


The jet mixing model 


Figure 


4-2. The air entrainment zone (e) 1s composed of a fuel rich 


zone Ce) 


where soot formation occurs, a products zone (d) 


where intense reaction occurs, and an outer zone composed 


primarily 


Mass meh) 


Sit cal Ty Outside the entrainment zone is the air 


(a). For the purpose of soot formation 


Calculations, a local equivalence ratio of unburnt fuel 2 


the fuel rich zone, Ru 1s defined. A series of equations 
ar@ proposed Oy Khan at al for @, calculations the 
significant points of which is that $F is inversely 
proportional Eo alr entrainment which 1s in turn a 


progressive function of air swirl 


1s not, 


mechanism of soot formation under high temperature and 


reaction 


combustion 


equation: 


ds 


al 


where: ds 
dt 
Vio 


aaa 


Lt 


Ama 4 
Ul 


time conditions characteristic of diffusion 


may be characterized by an Arrhenius 





soot formation rate 


volume of soot formation zone (c) 

cylinder volume 

partial pressure of unburnt fuel in soot 
formation zone (c) 

local temperature of soot formation zone (Cc) 
soot formation rate coefficient 

soot combustion activation energy 

universal gas constant 
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ratio and engine speed. ¢ 


Lt 


however, directly coupled to $6 charge. The overall 


short 


flame 


type 
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FIGURE 4-2: Soot production fuel jet mixing model of 
Khan et al [18] 
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After formation, soot particles coagulate and/or combust from 
a typical initial diameter of 250 A where combustion of soot 
coagulate 1s again controlled By an Arrhenius type equatioan 
meeportional to the partial pressure of O.- 

The values of constants and coefficients are, of 
course, a function of individual engine design and fuel 
structure selected to Gest "Curve fit" experimental data. 
Valuable qualitative conclusions, however, stand out from 
the work by Ehan et al: 

ve The rate of soot formation depends on three 

principle parameters: local F ee ee 
T, and 6 in the fuel rich soot formation 
zone. The effect of varying overall = 4$ 
charge 1s moderated as Di is always near 
unity. Swirl ratioa, however , has a 
significant effect. A comparison of soot 
concentrations measured by KFhan et al, 
Flynn et al, and Dent et al is shown In 
Figure 4-3 which graphically illustrates 
soot concentration, thus luminous flame 
radiation, suppression in high swirl Cl 
engines. 

a Approximately 60%, as measured By Khan et al 
Ci8], of soot initially formed is burnt as 
it moves from the cooler fuel rich zone (c) 
into the products zone (d) where soot 


oxidation incandescent radiant emission 1s 





a grey body, optically thick. The exhausted 
soot is that portion which escapes the high 
temperature products reaction zone (d) in 
the process of turbulent diffusion from the 
fuel rich to the lean zone. 

ee High speed photograohy indicates that the 

majority of soot formation occurs at the 
beqinning of the main combustion phase, 
within an extremely short period Csee 
Figure 4-3]. The total amount of soot 
remains constant in the later combustion 
stages. 

A clear conclusion of this model is that both soot 
formation rate and soot combustion rate are progressive 
functions of fuel and O., partial pressures. The New Cycle 
Engine, operating at significantly reduced expansion cylinder 
pressures, should therefore experience a marked reduction in 
luminous flame radiant heat transfer. Due to lack of (built) 
engine test data however, the degree to which O-. decreases 
relative to magnitudes observed in high ie CI engines must be 
estimated from indirect empirical relationships. The quality 
Of this estimation, based principly on reduced pressure 
consideration in low swirl conditions, must be tempered with 
the knowledge that many other operating parameters have a 


secondary etfect on flame emissivity and temperature. 
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FIGURE 4-3: Comparison of in~cylinder soot conentration, 
A» versus crank angle as observed by various 
investigators [10 
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ie SECONDARY FACTORS AFFECTING COMBUSTION SOOT 
CONCENTRATION 
Extensive tests performed by Flynn et al on aA four 
stroke, Girect injection CI engine, with piston Cup 
combustion chamber enlarged to suppress swirl ratio and 
enhance soot production, illuminate many of these secondary 
effects; the important ones of which are listed below (Cidi: 
Le Ge increased with inlet manifold pressure 
increase at nearly proportional rate; 
much less ignition delay and fuel/air 
pre-mixing occurs causing increased 
tendency for soot formation in diffusion 
flame. 
Ze AS injection time was advanced an 
increase im Tq and an increase in peak 
apparent radiant temperature occurred; 
due to fourth power effect a significant 
increase in peak radiant emissive power 
results. 
aoe Decreasing cetane number with fixed 
reference fuel families encouraged 
ignition delay, increased fuel/air pre- 
mixing, amd reduced soot formation; at 
supercharged operation this effect is 
effectively suppressed due to predominant 
effect of manifold pressure on ignition 


delay. 





4. Increasing C/H mass ratio fuels at fixed 


Gaetane nunbers encaduraged Ur, due to 
simple availability of carbon 
concentrations the rate is not 
proportional to C/H mass ratio and 


apparent only at high ratio aromatic 

fuels. 
Flynn et al additionally note that the apparent flame optical 
thickness, Keeley diminishes approximately in inverse 
proportion to increasing cylinder volume. This conforms with 
the Ehan et al model and Figure 4-3, indicating that the 
number of incandescent radiating particles in later 
Pemetotion stages 15 essentially constant. 
a 4 CHARACTERISTIC HEAT TRANSFER EQUATIONS 

Emission as a function of wavelength for clouds of 


small particles obeys a variation described by equation (55): 


ene 275 =n 
BE, = 1 - exp(—k,1p/a* 7”) (53) 


The spectrum is near grey, and for simplicity of calculation 


can be converted to equivalent grey body emissivity: 


€ 


I. 44,04 


n 
[ Gq, 5A 


where: an is range of significant A band 


The effect of combustion cylinder wall reflectivity, § as 


i? 


demonstrated by Oguri et al in Figure 4-1, may be corrected 


30) 





for by determination of an apparent flame emissivity, € 


according to equation (S35): 


where: Kf = apparent flame absorption coefficient 


per unit length 
When flame envelopes the entire cylinder volume, observed by 
Alcock et al Citd to occur approximately 190 to 20° after 
ignition start, I. may be calculated by eguation (47) as a 
function of instantaneous cylinder V and Aus Sitkei eat al 
CiSjJ state that surfaces of pre-swirl chambers, piston crown 
and cylinder heads have §) = 0.1 (approx.)3; cylinder sleeve 
Sy = 0.2 (Capprox.). It should be noted that estimates for %, 
vary greatly, in part due to different modeling of cylinder 
surfaces as specular reflective to semi-rough diffusive 


reflective characteristic of oil Soaked, soot-contaminated 


surfaces. Or, may then be determined by equation (34): 


— : 4 = 4 eases > 
Ge. = © or Au els ie J¢1 hre/26ood sec) (35) 
FeTuU/sec 
ca SITKEIT ET AL METHOD OF AFFARENT FLAME EMISSIVITY, €.1, 
CALCULATION 


A useful empirical relationship used by Annand C3] and 
Dent et al CiO] based on known soot concentration is stated 
Poeeqtkation (37) in U.S. units: 
ke = 5898 8 ft7! (57) 


where: 8 = soot concentrations in lbm/ft™ 


Si 





For estimation of © without knowledge of soot 
concentration, Slitkel et al propose the following definition 
of oS as a function of instantaneous F, $6, le and crank angle 


past ignition start C1isi: 


2 045K OF le, 


ean o Pia 4 (38) 


where: €, = emissivity of infinitely thick radiating 
flame; = 0.8 


ke = apparent flame absorption coefficient per unit 
length-unit pressure (atm) 


Values of EG are determined from Figure 4-4, at applicable 
crank angle, and excess air factor equivalent to oe Sitkel 
et al based Figure 4-4 on measurements conducted by a photo- 
Gelectric pyrometer unit through a quartz view window insert 
to a swirl pre-combustion chamber of a four stroke Steyr 
diesel at r. equal to 18.5. Recognizing that e¢_, is dependant 
on many additional variations i1n engine operating parameters, 
Sitkei states that to a first approximation [14, pg. S46]: 

“in the absence of other data, it is felt that & 

Given here is universally adaptable to al 

engines of the diesel type.... with the known F-V 


diagram of the engine under test the variation of 
Kk can be found" 


p 
Inherent in this assumption is an assumed %) of 19% to 20%, 
an ignition start near TDC, a rather arbitrary value of eG = 
O.8 based on observed steady state diesel fuel flame 
emlissivities in combustors, and the belief that iS observed 


in a 2.6 Cm pre-combustion chamber may be extrapolated to 


main cylinder emissivities. 








— 20 40 60 80 


+0 from ignition start : “degrees 


FIGURE 4-4: Variation of ky versus angular displacement for 
_ calculation; Sitkei et al model [14] 
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Tt should be noted that the values of ©. determined by 
the above method are generally less than observed By a number 


of other investigators in direct injection CI engine testing. 


Dent et al state (C19, p. 1760): 
"The essential results of Sitkei and Ramanaiah 
are in agreement with those of Flynn et al in 
terms of the trend of flame emissivity and 
optical thickness with crank angle, load, and 
overall air fuel ratios. However, the magnitude 
of the flame emissivity and optical thickness are 


approximately half." 


For first approximation of predicted peak &€_ in the New Cycle 


a 
Engine expansion cylinder, it is therefore believed that the 
&. obtained by the Sitkei et al method should be multiplied 
by a factor of Zz. Tt is also noted that Sitkei et al agree 
with typical Oe. / range of 20% to 20% even at reduced a 
Values due to Seine of Belinskiy’s significantly hiqher 
le profile shown in Figure 4-3, Vice direct optical 
measurement. 
ie 6 LOW COMBUSTION FRESSURE FOWER UNIT METHOD OF 

APPARENT FLAME EMISSIVITY, €., CALCULATION 

For first approximation purposes, an alternate method 
of predicting oo is to examine the characteristic emissivity 
equations used in power units operating at the low peak 
combustion pressures of the New Cycle Engine. Different 
properties, continuous flaming and non-reciprocating 
mechanics of gas turbines do not allow Sroad application of 
flame emission data to internal combustion engines. However, 


actual combustion kinematics and flame luminosity observed in 


power units at pressure conditions more reflective of the New 
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FIGURE 4-5: Measured Tp in CI engines versus crank angle as 
observed dy various investigators [10] 





Cycle Engine should not be ignored. It 1s therefore 
considered that a fair prediction of New Cycle Engine peak e . 
only can be ascertained by use of gas turbine power unit 
emissivity equations using diesel fuel properties. 

For gas turbine application, €_ can be determined by 


cA 


equation (39) C194]: 


$= i-e* CSF) 


where: z = -38636.55 Fl(fuel/air mass ratio)1,]1°9Tg !*%z 


The flame luminosity, #&, can be determined as a function of 
fuel C/H mass ratio from Figure 4-6. For commercial liquid 
fuels of conventional character ranging from gasoline to fuel 
Oil, C/H mass ratio can be estimated using equation (60): 


Eemass ratio = .16015p —- 1.46 (60) 
H 


where: p fuel density in lbm/ft” is referenced to 
Ss°R (288°) 


is 7 ROLTZMANN NUMBER FUNCTION METHOD OF FLAME 

Temperature, Tye, CALCULATION 

A consensus opinion of investigators support that for 
first approximation purposes, radiant temperature is equal to 
Tei various Tg profiles published from high rr CI engine 
testing are shown in Figure 4-5. An empirical correlation 
allowing calculation of Ny is the Boltzmann number function, 
Tg/T, = f (Bo), established in Figure 4-7. The Ho” is 


defined in accordance with equation (61) Li4d: 
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FIGURE 4-6: Flame luminosity,&, as a function of fuel C/H 
mass ratio in low combustion pressure power 
Paice oO. 
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Pal 


3 fuel 
Bo®* = - 7 (61) 
ae oe ae 
where: B = fuel consumption rate 
Cy fuel = mean total specific heat of combustion 


products of 1 1bm of fuel 
For ease of calculation, Bo and Cy fuel may be calculated for 


the New Cycle Engine using the following identities: 


(SFC) (Aneat) (Menarge? (Nrpm? (69 min/hr) 
B Se Coe) 
lbm/hAr fuel 2944 BTU/hr -HF 
are fuel x= (Air/Fuel mass ratio) (C, combustion products) ror) 
substitution of the above identities into equation (61) 
yields: 
(54) 
; C (SFC) (Oheat) (Me parge? (Neem? (Air/Fuel Mass COO ee 
Irn rnnn nnn 
aioe 


where: C = 1.2768x107 


Calculations of Ty, at a known Tg and €, from the Bo function 
requires an fe eb dwe process matching an assumed Vs both 
directly from the generated pea aie fraction and its coupled 
Bo* from Figure 4-7. An increasing of Will decrease the Bo , 
thus decrease [eee resulting in an increased calculated I[_. 
It is noted, however, that increases in Cs and le are not 
proportional; over a wide range of Bo* ( eli TQ/T 5 slope 
is only slightly positive. At Bo* greater than 1.2, 


characteristic of low luminous flames, the value of Ty, is 


essentially independent of Bo* (thus €,) and equal to Tg/0.8. 


oH 
Co 
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EG PREDICTION OF FEAK COMBUSTION €._ and Ty 


For calculation of peak ©. and Te in the New Cycle 
Engine expansion cylinder at ‘full load, the +foallowing 
ASSLUMPtTLiONs are made: 

ie The expansion cylinder reflectivity, § = 0.1 

es The flame growth encompasses the entire 

cylinder volume at 20° maximum crank angle from 


Ignition start Ccrank angle = 65°]: flame 


growth rate 1s constant. 


ss The flame luminosity, &@, calculation is based 
on medium speed diesel fuel C/H mass ratio = 6.3 

4. The fuel/air mass ratio = .0905 (1/20) for @ = .7%5 

Je Te at fuel C/H mass ratio = 6.3 may be 


approximated by equation (322) as a function of Tg 

felts CHa. fuels (C/H mass ratio = 6.0). 
The Sitkel et al method of €. calculation was developed on a 
test engine exhibiting typical high r. combustion pressure 
profiles, and TDC ignition start, both conditions not met in 
New Cycle Engine operation. The peak Ky from Figure 4-4 at 
the correct excess air factor value will therefore be applied 
at the crank angle of maximum combustion pressure where flame 
growth reaches 100% cylinder volume. 


{ 


Sitkel Method; C peak 








Qe 50 
“ta peak oo 
excess air factor = (.75)7+ = 1.33 
peak KG at 1.23 excess air factor = 0.63 


(see Figure 4-4) 


FO 





sing equations (47) and (38), plus applying a oS correction 


metecai of =: 


9p y Ar i ky 304K Plp f, 2, 
ate +t ¢t? ft (#teata) 7! 
65° 3.76064 488331 3.944 4210 63 304 210.420 


p = .~420 
eal os Ai 


at @ = 65° 


Low Combustion Fressure Fower Unit Method, ©. peak 


flame luminosity, &© = 4.0 
(see Figure 4-6) 


Using Equations (47) and (359): 


Q p To y A flawe growth = lork l C. 
she tR >t? 1 tt 

50 4.24684 2259.6  .309053 3,358 25 0782 «6382 318 

55 4.09172 2571.21 .366348 3,545 50 1757 455.366 

60 3.92401 2864.13 426316 3,741 75 2906 9.477.379 

65 3.76064 3136.02 488331 3.944 100 4210 480 =. 381 

70 3.26107 3064.15 551767 4.152 100 4518 447.380 


peak C . = ,361 


at 9 = 65° 
A repeated calculation for both methods at peak combustion 
oressure greater than SO atm, typical for direct injection, 
high ro Cl engines, results in peak €_, = 1.0, which closely 


agrees with the expected value (20.9). 


Roltzmann Number Function Method, peak ee 
c= O 


The known values required for calculation of the Ho* at 45 


crank angle, peak €, of (.420 + ~281)/2 = .40 are: 


yal 





Ps Tg = 3126.0 °R 

ee SFC = .284298 lbm fuel /HF—-hr 
: Meat = 300 BTU/IbMm charge 
Teharge = 024537 Lom 


we Newey = 612.3 


> ¢, = . 247+. 243x107 (S126.0°R) = . 33352 BTU 
S15 lbm—-°R 


ys air/fuel mass ratio = 20 
Sie Aut = 3.944 #t*~ 


substituting the above values into equation (54) yields: 


(1.3768x107) (.28429) (S00) (. 024537) (818.3) (20) (. S252) 


Bo* = = 
6.435x101% 
Cel 
From equation (33), an estimated peak o Of O.40 at an 
assumed §) = G.1 yields a peak €, equal to 0.426. The 


predicted peak Te may therefore be calculated from iteration 


+ 


using Figure 4-7 at Ho” equal to: 


6.425x%1019 = 1.511x10!! 


4267 5 ~ TT,” 


Assuming a value of Ve = Tg/O.8 = 3S920°R, the calculated Bo* 


equal to 2.1. This Bo* is greater than 1.2, therefore the 
estimated peak Ve equals i2920°R. Rounding off for first 
approximation: 


peak Vs = 4000°R 





To check the accuracy of this procedure, a peak © 
Value of 1.9 15 assumed commensurate with observed peak e in 


fPeeect injection, high a Cl engines: 


tO 10 


6.425x 10 = 6.4359x 10 


Py 
2 
° * 


f He 


te ee £ 


From Figure 4-3, the observed peak Ty, is typically equal to 
4=200°RS Tg/T now equals 0.735. Substituting 4200°R into Bo* 
= es. 10+ 7. yields a value of 0.87 which, entering 
Figure 4-7, extrapolates to ey ae = 0.76. The iteration 
essentially matches, indicating that a peak © and peak ee 


equal to 1.9 and 4200°R respectively, are in fact reasonable 


Values for direct injection, high r_ Cl engines. 





4. CALCULATION OF ORDER QF MAGNITUDE REDUCTION 

IN COMBUSTION RADIANT HEAT FLUX 

Woschni did not separately measure radiant heat flux in 
his engine testing; for ease of calculation he simply 
absorbed the radiant component into the forced convection 
heat transfer model by use of time averaged, “curve fitting" 
multiplicity coefficients. Woschni in fact believed that 
ieee p. ©0974 — 20751: 

"the total radiation amounts to no more than 20% 

of the total combustion term...In steady flames 

there are usually large combustion chambers with 

great volumes of flames, and because oar the 

small velocities there are small convective heat 

transfer rates only, so that in these cases 

radiation is responsible for an appreciable 

amount of the total heat. Because of the high 

velocity and density of gas, the convective heat 

transfer coefficients in internal combustion 
engines are about ten times higher, so that 
despite the higher temperatures the radiation 

Plays an appreciably smaller part." 

While the thrust of the statement is correct, Woschni 
unfortunately underestimated Qe, basing cursory radiant heat 
transfer calculations on the assumption of cL = 0.6 fram 
steady state flame observations and employing for the unknown 
radiant temperature, Tg, vice T,. 

For the purpose of this calculation, establishing 
Woshni peak &. and Ty values of 1.0 and 4200°R better 
conforms to published data on CI engine testing of similar 
design, direct injection, high ice characteristics. [It is 
noted that actual observed peak Gc, values are typically equal 


mem. 7 Vice 1.0. To maintain continuity with C prediction 


methods used in section 2.8, however, a value of unity is 
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assumed aS predicted for high peak pressure conditions. 
Without knowledge of the radiant component of dy and Gnce 3 
a function of crank angle, it is not practical to attempt 
calculation of a mean order of maqnitude reduction ratioa, 
Gyy/Iace> over the effective 990° crank angle domain where 
luminous radiant heat transfer can be significant Co == asymm) 
Pine gh ros medium speed CI engines. Fortunately, ain this 
case it is not required. Making the reasonable assumption of 
equivalent €., Ty, decay rates, flame emissivity in the New 
Cycle Engine expansion cylinder should certainly decay below 
O.1 at advancing crank angles from ignition start prior to 


Ba) Thus, peak Qy/dyce should be a minimum, conservative 


estimation of the mean. Therefore, from equation (3546): 


= Boe 
Gy Ow Celie Tw dw 
— CK peak « 7. 20 
INCE ONCE Galle “TW Once 
peak 
madi ation ie = GOOR 


Substituting in the respective CK and ee Values previously 


determined: 


Bi 1.0¢4200% 


minimum = 7 
INCE O,.4(49000 °-900) 


-9004) 


-9O47 approx. 
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CHAPTER Pave 


CALCULATION OF CORRECTED HEAT TRANSFER LOSS/HEAT 


Fpte Ase tmae 7 cl amine 


Total heat transter 10ss5, Ge, for a natural 
aspirating, four stroke CI engine includes heat transfer 
during the induction and compression strokes which do not 
exist ain the two stroke New Cycle Engine expansion cylinder 
cycle. In the induction stroke, the heat transfer driving 
potential, Ta-Tw, is negatives induction charge heating is 
boosted By residual gas mixing. In the compression stroke, 
Mifebatic compression may boost Tg to greater than 1500°R, 
From Figure 2-1, it 1s observed that a significant positive 
Tg-Tw occurs only in the last 30° of the compression 
stroke. The order of magnitude of this positive driving 
potential compared to that observed in the combustion- 
expansion and exhaust strokes 15s still quite small. It may 
therefore be assumed that the induction and compression 
stroke contribution to OQ for a natural aspirating, four 
stroke, high ie CI engine is not significant. A direct 
coupling of observed O.,/Q+ ratios may then be applied to the 
New Cycle Engine expansion cylinder. 

From Chapters 2, 3, and 4, the following principal 


Values have been determined: 


1. Ry uncorrected = 19.751%3 Woschni correlations 


“" 


a Ge /Q.- Woschni engine = WY; .2 = Y= 3 





3. Forced convection quy/dyce = 6-087 

4. Radiation qy/qyce = 3.047 minimum 

R+- must be corrected tor the failure of Woschni’s 
correlation, based solely on a theoretical forced convective 
heat transfer model, to properly scale down luminous 


radiation heat transfer to low ae Operating conditions. 


Correction to R+ = Rep + Re 
At high Ir operation: 


Che , OR / (Oheat) (Me parge’ : Re. ae 
a. (ile, +t: Geo? /fahest) (Me parge’ Re. a Rep 
Re, = YR. + Reo 
Re (i-W) = WRep - 
Re = Wo Rep (65) 


Ley 
Substituting into equation (65) the known order of magnitude 
reduction ratios for low re operating conditions and 


establishing equality: 


At low r. operations 


= are % Beas ¥% 
<n? @ O47 Re o— WY 6. O87 Rec 
1-W 
% Bee nto * 
Maxi MUM 1-W 
AS Woschni ’s correlation properly scales down forced 
convection heat transfer only, ee = (1 - W) T+ Therefore a 


conservative, first approximation prediction of correctad 
total heat transfer loss/heat release per cycle, ae 


\ 


for .2 2 We ais : 
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Forced convection heat transfer loss, ft 


a 
ot a ae _ = e 
iene = on = 15. 809% 


Radiant heat transfer loss, ie 


* 


= 27 OC aa 283 
Re = 2.00¢ WY Rep Rg GD) aS Ba Je W) it. 


1-W 1-W 
= 2.00 WY Ro 


= 2.00 WY 19.7617 


ye 7O4y7, = Te. = 11.8572 
Max l mum 


Total heat transfer loss, r-* 


H * a ea 


= (2. OOW+1—-W) RL 
= C1 +H) Rr. 
= (1+W) 19.7612 


ee Roe Pee oss 7. 


maximum 





CHAF TER SIX 


CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE STUDY 


A semi-empirical correlation for mean instantaneous 
heat transfer coefficient, developed by Woschni an 
conventional four stroke, direct injection, quiescent CI 
engines, was modified to predict heat transfer loss in a 
computer-generated preliminary design engine cycle. The 
engine divides the four stroke functions between a 
compression cylinder and an expansion cylinder interconnected 
by a static regenerative heat exchanger. Due to regenerator 


boost heating of the air charge, required compression ratios 


a the New Cycle Engine expansion cylinder will be 
Significantly reduced, and increased thermal efficiency 
achieved. To verify applicability of Woschni’s correlation 


to expansion cylinder heat transfer, points in the cycle were 
examined comparing Woschni’s correlation to fully developed, 
turbulent flow, forced convection relationships for flat 
plate and pipe flow. The mean instantaneous gas velocity was 
estimated from the construction of the charge impulse 
transient existing into the expansion cylinder from the 
regenerator port. An order of magnitude reduction in forced 
convection and luminous radiant heat flux was calculated to 
investigate the ability of Woschni’s correlation to scale 
down correctly total heat transfer to low compression ratio 


operation. 





It was found that the impulse velocity transient was 
sufficiently suppressed to approximate quiescent conditions. 
Tt was found that laws of similarity incorporated into 
Woschni’s correlation appear to allow uncorrected application 
at low compression ratio operation for forced convection heat 
transfer only. The order of magnitude reduction in 
combustion-expansion stroke forced convection heat flux was 
found to be as much as twice the order of magnitude reduction 
in radiant heat flux. A correction was proposed for the 
application of Woschni’s correlation to achieve a conservative 
prediction of heat transfer loss in the New Cycle engine 
expansion cylinder cycle: 

Gl. = 1.25 Ge. (67) 

corrected uncorrected 


where: Q, = a de 


[ 
uncorrected i 


3 
feta. 460 P°? rc, Sp3°9 a °* Tg «9 
Cy = 6.18 exhaust stroke 
A Ras combustion-expansion stroke 
This study investigated heat transfer in engine 


operating conditions significantly below the range for which 
published correlations exist. Theoretical work and actual 
measurement of an engine prototype remain to be completed to 
prave the validity of Significant extrapolation of 
reciprocating engine heat transfer correlations. This study 
points out three specific areas requiring future study 


concerning reciprocating engine heat transfer prediction: 


LOO 





Quantitative prediction of flow pattern magnitude 
and distribution as a function of crank angle and 
location. 


Quantitative prediction of soot concentration, 


apparent flame emissivity, and radiant 
temperature aS a function of operational 
variables in low pressure diffusion flame 


combustion. 

Quantitative effect of appreciable fuel 
injection and ignition delay in the expansion 
stroke of low compression ratio reciprocating 
engines on incandescent particle concentration 


and luminous flame radiation. 


er 


= = <- _ 
= 
a ——— 
> 
= => = 
- ]—_- 





(3) 


(4) 


(5) 


7? 


(3) 


Rei ene NCES 


Carmichael, A. Douglas, Frofessor of Fower Engineering, 
Massachusetts Institute of Technology, September 1983. 
Rohsenhow, W. and Choi, H., Heat, Mass, and Momentum 


POO TTT) PORTIS OED 6 had ETT FON ties TAO SOIREE OF 19 CD ff 0 RAI 


Transfer, Frentice-Hall, Inc., Englewood Cliffs, N.J., 
PO). 

Annand, W., "Heat Transfer from Flames in Internal 
Combustion Engines", Heat Transfer in Flames, N. Afgran 
and J. Beer, eds., Scripta Book Company, Washington 
Petweto/s, Chapter 24, pp. 377-390. 
Hohenberg, G., "Advanced Approaches for Heat Transfer 
Calculation", S.A.E. Transactions, Vol. 88, 1979, Faper 
P7OosesS, Ppp. 2788-2806. 
LeFeuvre, 7T., Myers, F., and Uyehara, O. "Experimental 
Instantaneous Heat Fluxes in a Diesel Engine and Their 
Correlation", S.A.E. Tramsactions, Vol. 78, 197697, Faper 
690464, pp. 1717-1729. 

Annand, We, "Heat Transfer in the Cylinders of 

/ 


Reciprocating Internal Combustion Engines", EEoe . 


Mmiotmea Neem. Emgrs., Vol. 177, 1963, pp. 9732-979. 


Mansuri, S. and Heywood, J., "Correlations for the 
Viscosity and Frandtl Number of Hydrocarbon-Air 
Combustion Froducts", Combustion Science and 


fecamoledy, Vol. 25, 1780, pp. #xsl-250. 
Oguri, T. and Inaba, 5., "Radiant Heat Transfer in 
Diesel Engines", S.A.E. Transactions, Vol. 81, 1972, 


Paper 720023, pp. 127-142. 


ace 





(9) 


i > 


oil) 


Ci>) 


C13) 


(14) 


(15) 


Woschni, G., "A Universally Applicable Equation for the 
Instantaneous Heat Transfer Coefficient in the Internal 
Combustion Engine", S.A.E. Transactions, Vol. 74, 1967, 
mapem 670951, pp. 2065-3084. 

Dent, J. and Sulaiman, S., "Convective and Radiative 
Heat Transfer in a High Swirl Direct Injection Diesel 
EnGgime », O-Ae&e Transactions, Vol. 86, 1977, Faper 
7/7/0407, pp. 1758-1783. 

Alcock, J. and Scott, W., "Some More Light on Diesel 
Combustion", Froc. Instn. Mech. Engrs., Automotive 
Pav tsieom Froceedings, 1962-1963, pp. 179-200. 
Fitzgeorge, D., and Allison, J., "Alr Swirl in a Road 
Vehicle Diesel Engine", Froc. instha | Mech. . Sagrms.. 
Automotive Division Froceedings, 1962-175. oe. LSA 
7 Se 

Heinein, N., "Instantaneous Heat Transfer Rates = and 
Coefficients Between the Gas and Combustion Chamber of 
a Diesel Engine”, S.A.E. Transaction, Vol. 74, 1965, 
Faper 650096, pp. i-1il. 

Sitkel, Sax, Heat Transfer and Thermal Loading 
in Internal Combustion Engines, F. Szoke transl., 
Akademial Fiado, Budapest, 1974. 

Sitkei, G. and Ramanaiah, G., "A Rational Approach tor 
Calculation of Heat Transfer in Diesel Engines", 5S.A.E. 


Transactions, Vol. 81, 1972, Faper 720027, pp. 165-174. 


10.5 





(16) 


a7) 


eo) 


(19) 


Flynn, F., Mizusawa, M., Uyehara, O. and Myers, F., "An 
Experimental Determination oF the Instantaneous 
Fotential Radiant Heat Transfer Within an Operating 
Diesel PAT eis Sees  HeANSACClIOM, © VOL.  O1,0 °° 1972, 
Rage 720022, App. FI-126. 

Chapman, M., jFFriedman, MoM. and Aghan, A., "A Time- 
Dependent Spatial Model for Radiant Heat Transfer in 
Diesel Engines", S.A.E. Transaction, Vol. 92, 1982, 
Peper BS1l/725, po. 13-20. 

Khan, I. and Greeves, G., "A Method for Calculating the 
Formation and Combustion of Soot in Diesel Engines", N. 
Afgran and J. Beer, eds., Scripta Book Company, 
Washington D.C., Chapter 25, pp. 291-403. 

Eretschmer, D. and Odgers, J., "A Simple Method for the 
Frediction of Wall Temperatures in Gas Turbines", ASME, 


1978, Faper 78-GT-90, pp. 1-7. 


1O4 





AFFENDIX A 


The purpose of this appendix is to present preliminary 
design parameters and the computer-generated thermodynamic 
firing cycle for the New Cycle Engine expansion cylinder on 


which the heat transfer calculations are based Ci]: 


Freliminary Desiaqn Farameters: expansion cylinder 


i. Mean piston speed, Sy = 000 te7min (SS, S20teesec) 

Ze square configuration; stroke = bore diameter (L=E) 

SUE Neglect clearance volumes Vc = 9 

4. Neglect cylinder sleeve clearance 

ws Regenerator port exit cross-section area, Arn = 
0.80 €t 7" (1/3rd cylinder circumference x .75 inch 
depth) 

o Compression cylinder supercharqe pressure ratio=1; 


matural aspirating 


7 Bulescent fuel injection at nominal <£5° crank 
angle 

S. Taqnition start at nominal 45° crank angle 

ie Specific heat release per cycle, Ahneat = FOO 


BTU/1lbm charge. 


19. Specific fuel consumption, SFC = .284298 Llbm/HF hr 


fi. Fuel base = nominal CLHa, medium speed diesel fuel 
i@. Full load equivalence ratio, 6 = ./5 
13. Full load air/fuel mass ratio = £0 


14. Mean heat transfer surface wall temperature, 
Tw = 9O0°R 


1S. Specific gas constant for air charge and 
combustion products, R = 32.90 Llbf-fFt/1bm—-PR 


16. Air charge Frandtl number, Fr = .70 
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TABLE A-1 


New Cycle Engine expansion cylinder computer-generated 


thermodynamic firing cycle model as a function of crank 
angle; S° increments C1] 
crank angle Pressure mean gas temp. Volume 
2 Fy, atm tien V tae 


EXFANSION-INDUCTION CONDITIONS 


oJ Beare e & ~149.04 oe 4067E-O38 
10 eee) | 2137.29 of) eres. 
tS 2291593 Zone aia 7 ees = As. 
20) 4.82479 eee a ee AS 
ad 3267408 ny ae. k ~9821858 
oO Sa Lol i 2194.29 ~11849:3 
aya geooe lS Sever. t al ioe 
4) Sees & OE) ores 7 ~2£04916 
45 4.142883 LeoG.e . 2045 
EXFANSION-COMBUSTION CONDITIONS 
50 4.24584 2 ee. 22099053 
a ee ley ae eet 2 266548 
60 S752 40 | =~864.13 »~426216 
6g) oe 7O054 Seo. Oe »~488251 
70 See) 7 2064.15 oe eae. 
Ts =.288976 S99 7S. S&S 616005 
B80 Sn Sores alt ~. 680446 
8x ee DOS 4 2-984.735 »~ 744518 
SO ee 7 OSG Zee 1S feo 7 ee 
oS 1.95628 ess - 869441 
100 1.81064 2748.78 a FE9S42 
105 1.68741 Eee . 78698 
LG flo? 7 26779 1.042 
i a ee 1.492461 2646. 0358 1.99408 
120 1.417435 Lona! 1 to S 
125 pene o> So ee a 1.18847 
1:20 ~2F661 2567. .= 1 SOS 
Lev 1.24904 =~JI48.44 1.268357 
140 1.208357 2227 eo. Leos ot 
145 fe eo oer eect 
150 1.145635 2498.51 eer 
foo Lee Lo 7 2483.62 1 ee ie ae 
160 1.10286 =474.46 1.40061 
165 1.08796 ~464.96 1.41492 
170 tO 7 7 Oe 2457.07 lore oe 
173 1.06986 2450.74 1.43129 
BDEe 1.905626 2445.95 ieee aa 
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EXHAUST CONDITIONS 


RDC 
ros 
190 
eo 
200 
Loa 
210 
215 
mel) 
fae ct 
aol) 
Faia 
240 
245 
aot) 
eae) cs 
26) 
Oe 
270 
pes 
=o 
es 
2FO 
200 
2.0 
o2lQ 
Rei ad 
Sc 
a 
2.0 
a 
240 
24S 
ot) 
cra aad 


TDC 


P 


TABLE A-1 


1.035 
fac) 
1.035 
1.05 
ete Cs 
een) 
Heer Gr 
1.05 
eechss 
1.035 
1.035 
1.03 
anole 
1.035 
era 
1.05 
ett 
1.03 
eae 
iss 
1.05 
Lo 
1.05 
ier es 
ee Cpa 
1.035 
1.05 
eee) 
1.035 
1.058 
eee! 
1.035 
1s oa 
1.05 
1.03 
1.03 


atm 


107 


(CONT ‘D) 


Ge iN 





2459.06 
2434.65 
2420.7 

2426.81 
2422.99 
2419.13 
2415.26 
~411.64 
2407.99 
~404.41 
~400.9 

297 eae 
2294.14 
Si erty 
Soy aes 
2294.71 
EB Si 7. Fd 
227 Se a 
276. 24 
oT ees 
2a7 Pais 
Zoos s os 
ey Olerahes 
2564.46 
ene e 44 
22O0. 4 
eS Seay 
eo? «Oo 
foudde G1 
Roane — Ste 
Soe 41 
cis eae) A 
249.67 
2248.36 
2547.09 
2045.85 


244,55 


Veet 


te oe Se 
Wes La 
1.42515 
1.41492 
1.40061 
hoe og oe 


1.25978 


e . ~ ° : = .. 


@ an? cont tant oy? 


1.20292 
Peas 
Bee cea 
1.i8849 
Leas 
1.0941 
1.04201 
5 ag) et 
eg ee 
~8624358 
ee 
© 744526 
690464 
»~616023 
Sod 7es 
488248 
See eye coe 
© 206.564 
© 2O9069 
eae eis, 
Lf O4G29 
~1592238 
~118504 
Beesley st 
5 Oey 
JOS 04255 
ON aa 


oe 4O8SE—-O8 


() 





Pak ATM 





TOC 40 80 20 l6O 
CO s “degr ees 


foe reulomeor expansion cylinder firine cycle 


eeeGURE A-1; 
motoring P, 
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APPENDIX B 


The purpose of this Appendix is to present the 
derivation of several geometric relations for ae simple 
geometry reciprocating cylinder. Figure H-1 illustrates a 


standard cylinder with defined piston linkage symbology. 


1. Instantaneous piston speed, Sp 


by the law of cosines: 


om 
limes 


Z 2 P 
l = + 4 + 2ascose 


— = s 

s“ - (Zacos@)s + (a“-1_“) = 0 

by the quadratic formula, solving for s: 
—-(-Zacose) = [ (Zacosa)< = pani h i = 


peg z= Bs ae 7S 
acoss + Cl, + a~ (cos7~86-1)] 


S 


substituting in the trigonometric identity, cos*@-1=sin“9, 
and taking the positive real root: 


~™ ™ 


Ss = acoss + ch = a*sin*e1!/4 (Bai) 
differentiating s with respect to time, t: 
ds de ds 


Sp = — — . 
de dt de 


II 
T 
| 


“Nt 


. a~singcose 
Sp = 9C-asines - ———————————____] 
(1 *-a“sin*e) 1/4 
a 
SeuOeerctuting im the crank ratio identity, R-_ = 1 fa, and 
Simplifying: 


7 cose 
(Rum-sinte) + < 


substituting in the 6 identity, 9 = 2TNecpy/O0: 


LOY 





ae TNe ny Sine cose 
sp = —————————— £— i + SURE aiiEEEEEEE ener (B-2) 
60 (Ru“-sinta) 17 


2. Mean Fiston Speed, Sp: 


= T 
oP [ Sp dea 


ll 
a|- 


ih a aS TNe ysl ne cose 
= [ a bl +}——--________. i 73 pia 
w J 60 (Ro “-sinte) 1/4 = 
—2aNeeny rt” < sinecose 
- pe ae eyes eoraee fe 
60 Je Ja (RL “-sin*e@) 
—2aNec cing oo Ssinegcose 
6Q bh (Rosin ©) a 
Let u = cose and substitute into the integrand the equivalent 


ldentities: 


duo = -sine 
sine = 1l-cos9 = i-u= 
cos(m) = -1l 
cos(O) = 1 
—2aN —S 
- REM —1 
So = ——— CC 2 + [ aT du J 
60 Jit (Ro#-14ue) 17 
—ZaNn 
es +1 
= ee ae, - + [ ~ a7? du J 
60 Ji, (R *-1+u*) 
—2aNn +1 
REM S oS a 
= ee ES a ee Sy 4y2yt/s ) 
e 
63 —1 


PO 





ce 


ae 


op = Caro) 
40) 

- NR EM 

S5 = ———- 
60 

substituting in the crank radius identity, a = L/2: 

7 2UNeem 

sp = ———— ee 
60 


2. Lnstantaneous Cylinder Volume, Ve 


Neglecting clearance volume, V and sleeve clearance: 


c 3 


i 
= 
7 

tJ 
- 


V (BR-4) 


4. Instantaneous Heat Transfer Area, Aut? 


Neglecting clearance volume, V and sleeve clearance: 


c 3 
4 
where: L. = 1] t+a- s 


S 
Ss from equation (B-1) 


SPeeealctilation of New Cycle Engine Expansion Cylinder 
Freliminaryv Design Farameters: 


Deeelog Ly ay, Ky Newco 9 max Sp 


From Appendix A, the following preliminary design conditions 
are known: 


ee 


a. Clearance volume, ve = QO 3: L., Sots, bce SiS 
b. Sleeve clearance = 0 : B = bore and piston diameter 
c. Square configuration : Bz L 


d. Sp = 335.33 ft/sec 
1.43323 ft> 


fe Voryo = 9.80770 a 





Vepo = 1.45333 ft" = wR“ L, = wE™ (B-4) 
4 4 
B=tL=C 4(1.4222%) 11/75 


T 


RB>EeL = 1.22204 Ft 


II 
rm 
“sy 
b) 


From the crank radius identity, a 
a = 9.61102 ft 


From equations (B-1), (B-4), and the geometric relationship, 
Ly, =1l,_-+a-s at a given crank angle @: 


a 
ay Zon cont A) nr 
—,=1, + a - Cacosé + (1_“-a*sin“e) '/=] 
TR* 


which at known V, &, and a values at 9 = 9O° simplifies to: 


4(.8077) . me a 
= 1+. 61102-£.61102c0890° + (1, *-.61102"sin=902) */#3 


~ 
-—— Sai S ~906024 = se — £3755 
I elosesl, + IN6O24 Pee 7 2 4BG& 


lL. = 2.4440 ft 


prom the crank ratio identity, RL = l,/a: 


: = Aa 
Re 4.0 


From equation (B-3), Neey = 60 Sp/2c: 


60 (20.55) 
N 


REM 


2\1l.22204) 


Neem = 818.3 


From equation (8-2), recognizing that maximum Sp occurs at 


@ = 90°, 270°: 


eB Pram ") 0 C)? 
: ae TN ec inFo = cos so 


Piz 





> = alntn..- 
P imax Soa 


Gaol (2m) (Sle. 3) 
60 


| dees fC/SeC 


[3 
y 
i 





de 





BORE DIAMETER 


~—— B ——> 
y | = i 

; 4____ ke GAS VOLUME HEIGHT 
4, 4 | =l.+da-S. 


: (Ver O 


1, CONNECTING ROD 


mw SS ee SS 


a CRANK RADIUS 





\ 
rs; 
/ 

yi RANK SHAF T : 
ee R eee 
ome 


CRANK RATIO 


Pega o—-is ooandard reciprocating cylinder and piston 
linkage geometry 
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AFFENDIX C 


The purpose of this appendix is to present formulas for 
calculation of essential thermal properties of principle 
gases asa function of mean gas temperature, Tq, and lean 
Hhydrocarbon-alr combustion products, CHa. fuels, as a 


function of Tg, pressure, F, and equivalence ratio, $@. 


Es maomeci ple gases 


TABLE C-1 


Coefficients for calculation of C_ and 4 of principle gases 


as a fumction of Tg(°R) [6] i 
Valid for Tq + 450°R 

Substance AMA YG, mara NONE OA ju. 
air ao 2224 eee Ley, - OO Zao 
nitrogen No 23 ee 2478 aren) iki Aa: 
carbon monoxide CO 3 eee 2a Oe ae oe aoa 
carbon dioxide CO. 44 ~141 eo ee —-2. 4867 1.704 
water vapor H.0 18 »~479 -,. 448% A389 = 1 1.499 
AMW = Atomic Mass Weight 

Me ease co 

lbm/ft-sec 
a 7 a 2 =o 
Cc, = Yo, + Y,I9g + YalTg? (G2) 


BTU/lLbm—-?R 


~~ 


2. Lean Hydrocarbon-Air Combustion products; C le fuels C7] 


wee tor YOO" R = fg = 7200°R and @ = 1.0: 


1.470x107 7 (Tg)* 
= J ee —) 


Me 
lbm/ft—-sec 1+.0276 





(C-—4) 


O.05 + 4.2°0y~-1) - 6.709-104 


Fir = 
¥(Tg) = Figure C-2 
Valid for 900°R £ Tg = T600°R, F = 1 atm, and 6 = .75 
4g (=e) 


Naa 


(see Figure C2) 
oat7 + 


Up 
ETU/lLbm-°R 











Saute 
e—e—e {00 





500s 1000 1500 2000 2500 
(900) (1800)  — (2700) | ca) (3600) (4500) 


Gq K ( 


PIGURE C-1: Ratio of specific heats, ie Or lean Sie wee leone 
Mioomeusuaon products as a function of FP, 12, 
and @; fuel C,H»5, [7] 


el 
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» 
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— 20 
I 3< (478) 
oo 
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< 

fo 
oO 10 
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[meune Ce2: ©. of lean hydrocarvon-air combustion products 


as a function of P, Tg, and Ms fuel CyHo, [7] 
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Thesis 
J8245 Jorgensen 


Predicted heat trans-. 
fer loss in the expan- 
sion cylinder of a two 
cylinder cycle reciproe 
cating engine. 


Thesis 
ie245 Jorgensen 
Predicted heat trans— 


fer loss in the expan: 
sion cylinder of a two 
cylinder cycle Peeipro— 
cating engine. 











